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HASK 1; SPECTRQPHOTOMETRIC OBSERVATIONS AND 
ANALYSIS OF THE ASTEROIDS 

(Principal Investigator; Clark R. Chapman) 

General Remarks 

Daring the final quarter , progress was made In several areas on 
the apteroid project. The report must he lirief because the P.I. is 
travelling out of the country . 

Observing Program 

The last observations were conducted at Lowell Observatory in 
late May aiid were reported in the previous report. The next observing 
runs will be in mid- and late-September at KPNp, on the 1.3 and 2,1m 
telescopes, respectively . A selection of targets was assembled, con- 
centrating on the important faint objects that will be available — the 
unique asteroid/comet Hidalgo and 1620 Geographos. 

Reductions of Previous Data 

A trip by Chapman to M.I.T, resulted in substantial progress 
towards a final reduction of data taken during observing runs over the 
past year and one half. A calibration for ee Lyr through the moon to the 
sun was settled upon as well as relative colors between a Lyr and 
various standard stars. The one remaiaing problem is to run appro- 
priate interpolations to the wavelengths of the filter-detector combina- 
tions actually used . 

Comparisons between Gaffey's and Chapman’s reduction programs 
showed that they were yielding the same answers. Procedures to 
corx'ect some of the data for minor digital glitches and for small coin- 
cidence problems were agreed upon and processed data have now been 
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received by Chapman. In view of the fact that P.S.I. will be 
obtaining a new, upgraded computer facility by September that will 
enable us to handle data on tape (rather than only keyboard entry), 
it was decided to postpone the final processing until Chapman’s return 
from the I . A . U . Colloquium . 

TH3AD Cooperative Studies 

The Tucson Revised Index of Asteroid Data (TRIAD) is now fully 
operational, with most major data sets received. The first major effort 
was to design observational criteria to permit classification of asteroids 
into the broad G, S, M, etc. classes on the basis of minimal data (e.g. 
ITBV colors). This joint project between Bowell (Lowell Observatory), 
Gradie and Zellner (University of Arizona), Morrison (KPNO), and 
Chapman (Planetary Science Institut j has resulted in a sophisticated 
computer program that classifies any asteroid on the basis of the com- 
plete suite of observational data, A manuscript is partly written. 

Chapman also assisted a fundamental application of the TRIAD 
data bank: performing bias corrections and obtaining distributions 
of asteroids by size, semi-major axis, and so on. The preliminary 
results (Figure 1) differ in some important qualitative ways from the 
previous bias-corrections in Chapman, Morrison, and Zellner (1975). 

The C asteroids evidently do not increase in frequency toward smaller 
diameters as rapidly as had been expected- Consideration of the Palomar 
Leiden survey (PLS) distributions at still smaller diameters suggests 
that the S objects may have a steeper fragmentation ’’tail” than do the 
C objects, although that hypothesis is not yet established . 

Lyon Meeting 

The I. A. U. Colloquium on interrelations among asteroids, comets, 
and meteorites is being held August 16-20 in Lyon, France . Chapman 
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-coneiuded preparations for this meeting at isrfaich he will present a 
TSTiew of the nature of asteroids as meteorite parent-bodies. 

Supplementary Color and Spectral Data 

Preliminary reductions of spectrophotometry, plus color esti- 
mates based on pre-1960 photometry, have been provided (see Table I) 
for those asteroids not already observed by Bowell’s on-going UBV 
observing program. 

Betulia Campaign 

The coordinated effort to observe Betulia in May was quite 
successful and succeeded in demonstrating that Betulia is a very unusual 
Earth- approaching object, perhaps a dead comet nucleus. Chapman 
prepared the following summary for publication in Sky and Telescope: 

On May 23rd, the tiny asteroid 1580 Betulia passed 
within 0. 13 A.U. of the Earth. Around the world, teams 
of astronomers attempted to follow it as it moved south- 
wards by more than six degrees a day. The results from 
this Betulia observing campaign will be reported in the 
journals during the next year, but already the preliminary 
data indicate that Betulia is an unusual Earth -approaching 
asteroid . It appears to be the best candidate yet for being 
the nucleus of an extinct comet . 

Ben Zellner and Jon Gradie of the University of 
Arizona measured the changing polarization of Betulia as 
a function of phase angle and from that determined the 
geometric albedo to be roughly 6% or less, making Betulia 
the blackest known asteroid inside the main asteroid belt . 

Also, at the Mt. Lemmon Observatory near Tucson, Glenn 
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Table I. Supplementary Color and Spectral Data 


20 spectrum is anomalous S 

31 spectrum is C; color is 
consistent with C or M 

36 spectrum is C 

38 color is intermediate 

57 spectrum is S; color is S 

64 spectrum is consistent 
with M or E 

70 color is consistent with 

C or M 

71 spectrum is S 

77 color is S 

92 color is intermediate, 
probably G 

110 spectrum is anomalous S 
(note inconsistency with 
UBV data) 

116 spectrum is S 

124 color is S 

153 spectrum is C 

154 spectrum is C 
179 color is S 

185 spectrum is C 
216 spectrum is M 
258 color is S 
306 spectrum is S 
344 spectrum is C 


345 color is consistent with 
Cor M 

372 color is consistent with 
C or M 

385 color is S 

386 spectrum is C 

415 spectrum is possibly M 
426 spectrum is C 
434 spectrum is M 
471 spectrum is O 
485 color is intermediate 
558 spectrum is M or C 
560 spectrum is C 
639 spectrum is S 

675 color is S 

676 spectrum is possibly C 

747 spectrum is C 

779 color is consistent with 
C or M 

785 spectrum is M 
801 spectrum is possibly C 
1036 spectrum is S 
1173 spectrum is C 
1580 spectrum is C 
1636 spectrum is anomalous S 


- 6 - 


Veeder and Marcia and Larry Lebofsigr measured the thermal 
radiation from Betulia at 10 microns. From their simul- 
taneous visual photometry, they hope to determine the 
albedo as well. 

UBV color measurements by Zellner and by Ed Tedesco 
at New Mexico State University' Observatory confirm that 
Betulia resembles the neutral-colored, low-albedo ’'carbon- 
aceous-type" asteroids — the first of nearly a dozen well- 
observed Mars- or Earth-crossing asteroids to do sn 
Narrowband filter photometry obtained with the 72 -inch 
Perkins telescope at Lowell Observatory by Clark Chapman, 
Ted Bowell, and Tom McCord, confirm that Betulia has a 
carbonaceous-type spectrum. 

Measurements of Betulia’s lightcurve were made by 
Tedesco, by Bob Millis at Lowell, and also by a team in 
Perth, Australia, coordinated by I. Hikoloff. Tedesco's 
reductions of some of these data yield a synodic rotation 
period on May 22 through 24 of 6.138 hours. The lightcurve 
amplitude was only 0.2 to 0.28 magnitudes, rather modest 
for a small asteroid, suggesting a reasonably spherical shape. 
The last result is consistent with radar echoes obtained on 
two nights at Arecibo by Gordon Petfcengill, Don Campbell, 
Brian Marsden, and Irwin Shapiro . 

Because of its low albedo, the polar imetric observers 
calculate that Betulia is larger than had been supposed — about 
6 km diameter . The strong radar signal and its doppler 
frequency spread, comhb.ed with the measured rotation rate, 
yield a lower limit to Betulia 's diameter of roughly 6-1/2 km. 
The radar data further reveal Betulia to be rough at the scale 
of the radar wavelength (12 centimeters) or larger. 
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Although nobody really knows what the nucleus of a 
dead comet is like, a carbonaceous composition seems 
cosmochemically more plausible than the minerals inferred 
to eslst on the other observed Earth-approaching asteroids. 
It is noteworthy that, of all the numbered asteroids, Betulia 
has the largest orbital inclination (52®^) and one of the largest 
eccentricities ( 0 . 49 ). Another asteroid, 944 Hidalgo, has 
an orbit even more like that of a short-period comet, ex- 
tending from the inner edge of the belt all the way out to 
Saturn; it will be the target of a new observing campaign this 
autu3xm. 
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TASK 2: LUNAR OCCULTATIQN 07 ASTEROIDS 

(Principal Investigator: Donald R. Davis) 

A search v?as made to find possible occultations of bright 
asteroids by the moon in 1976 and the first part of 1977, which is as 
far into the future as asteroid ephemerides have been published to 
date . All asteroids brighter than magnitude 12.5 at opposition were 
treated; there are 88 such asteroids within the next year. Of these 
88 candidates, 63 do not intersect the track of the moon in the sky* 
From the 25 asteroids which do intersect the lunar path, 17 are likely 
to result in lunar occultation visible someplace from the Earth's 
surface . Table II summarizes these potential events . These are 
preliminary predictions of the occultation events and more detailed 
predictions of locations and timing need to be made; David Dunham 
has, in the past, refined our preliminary predictions and announced 
the events to the astronomical community. A copy of the list of 
preliminary occultation predictions has been sent to Dr. Dunham of 
the Cincinnati Observatory . There are coii^uter programs now 
available in Tucson that allow detailed prediction of occultation events 
of this type. Although beyond the scope of the current task, future 
asteroid predictions should include some detailed predictions for 
specific events of interest . 
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Potential Occnltations of Asteroids by the Moon 


Asteroid 

Date 

Date and 
Brightness 

Lunar Phase 

Approximate 

Location 

4* Vesta 

7/24/76; 0^ 


3 days before 
new 

Mideast, Russia 

7 iris 

6/15/76; 16^ 

6/11/76; 10,4 

3 days after 
full 

Russia, China, 
Japan 

14 Irene 

6/11/76; 18^ 

6/11/76; 10,3 

1 day before 
full 

Russia, East 
Europe, India, 
China, Japan, 
Australia 

15 Ennomia 

7/14/76; 3^ 

7/11/76; 9.6 

1 day after 
full 

South Africa, 
South America 

22 Ealliope 

11/8/76; 20^ 

11/8/76; 10.8 

2 days before 
full 

Europe, Russia, 
South Africa, 
India 

59 Elpis 

8/U/76; 

8/27/76; 11.6 

1 day after full 

Australia 

63 Ausonia 

6/5/76; 18^ 6/1/76; 12,0 

Possible on 5/9/76; 3". 
Potential in July; no ephemeris. 

First quarter 

South Africa, 
Australia, 
Russia, East 
Europe 

80 Sappho 

11/4/76; 22^ 

10/16/76; 10,3 

2 days before 
full 

South Africa 

89 Julia 

10/1/76; 22^ 

9/29/76; 11.4 

1 day after 
first quarter 

England 

101 Helena 

9/9/76; 19^ 

9/20/76; 11.6 

1 day after full 

South Africa, 
Australia 

110 Lydia 

10/11/76; 24^ 

‘ 11/16/76; 11.8 

4 days after 
fuU 

South Africa 

192 Nausicaa 

6/5/76; 19^ 

6/1/76; 13.4 

First quarter 

Australia 


* Listed in A.E.N.A. , 1976 
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Asteroid 

Date 

Table II (Cont.) 

Date and 
Brightness 

Lunar Phase 

Approximate 

Location 

245 Vera 

12/4/76; 6^ 

11/12/76; 12.5 

2 days before 
full 

US, Western 
Hemisphere 

258 Tyche 

11/4/76; 11^ 

10/11/76; 11.4 

2 days before 
full 

West US, Hawaii, 
Australia 

409 Aspasia 

6/10/76; 4^ 

5/18/76; 11.2 

1 day before 
full 

US 

471 Papagena 

1/2/77; 22^ 

1/2/77; 10.9 

3 days before 
full 

Europe, England; 
Russia 

532 Herculiiia 

1/30/77; 16^ 

1/27/77; 11.9 

3 days after 
first quarter 

China, Japan, 
Russia 
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TASK 3: ANALYSIS OF SATXJHN*S KINGS 

(Principal Investigator: Michael Ji Price) 

Analysis of the visual photometric function in terms of anisotropic 
■scattering within the rings has been completed. Laboratory comparisons 
indicate that mineral contamination of the HgO ice particles cannot 
exceed 5 percent by weight- Results of the photometric study are con- 
tained in a paper, entitled "Anisotropic Optical Scattering Within 
Saturn's Rings, " recently submitted for publication in Icarus (Appendix A) 

Co-Investigator W, K. Hartmann has continued his experimental 
investigation bn rock and ice fragmentation. This material has been 
factored into a paper on Planetary Rings prepared in collaboration with 
Richard Greenberg (University of Arizona, Lunar and Planetary Labora- 
tory; transferring to PSI 8/16/76). The manuscript, entitled "Size ' 
Distribution of Particles ia Planetary Rings.," has recently been sub- 
niitted for publication in Icarus (Appendix B) . 
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TASE: 4: PROBING THE GIANT PLANETS WITH THE 

RAMAN EFFECT 

(Principal Investigator: Michael J, Price) 

Work on developing the Raman probe technique to investigate 
the physical structures of the outer planet atmospheres was completed. 
A paper entitled ”Gn Probing the Outer Planets with the Raman Effect" 
has been prepared (Appendix C) . It was recently submitted to 
Reviews of Geophysics and Space Physics for publication. 

Work on studying the limb-brightening phenomenon on Uranus, 
carried over from Contract NASW-2718, was completed. A paper en- 
titled "Limb-Brightening on Uranus: The Visible Spectrum" was sub- 
mitted for publication m Icarus (Appendix D) . It will appear in the 
1976 August issue. 
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“TASK 5: MARS DATA COLLECTION ANALYSTS . 

(Principal Investigator; William R. H^tmann) 

During the final quarter, the observations of Mars clouds, compiled 
^during several visits to Lowell Observatory in the last two years, were 
reviewed. A draft manuscript prepared earlier was extended and re- 
vised into the form of a paper to be submitted to Icarus for publication. 
'This paper describes and completes the project and is included as 
Appendix E. Early results were given in a verbal paper at the 1975 
meeting of the Division for Planetary Sciences, American Astronomical 
Society. 
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ABSTRACT 


Visual photometric function data for Saturn’s rings are analyzed 
m terms of elementary anisotropic scattering radiative transfer models 
'wMch involve the Henyey-Greenstein function. Limits are placed on 
the combinations of single scattering albedo and back -^scattering 
directivity which are permitted by observation. Particles with both 
microscopic and macroscopic lunar-like scattering properties are 
excluded by the analysis. Results are consistent with the ring particles 
being nearly pure spherical conglomerates of H^O frost. 


I. INTHODUCTION 

Interpretation of photometry of Saturn^s rings remains contro- 
Trersial (Follack, 1975). Specifically, the relative significance of 
■primary and multiple scattering is a matter of debate. So is the degree 
to vrhich mutual shado-wing occurs among the individual particles . 

Becently, Price (1973, 1974) rediscussed the optical scattering proper- 
ties of Saturn’s rings. E-vidence indicating that the photometric function 
at visual wavelengths is dominated by primary scattering, and that 
mutual shadowing is an irrelevant concept, was first presented. A-vailable 
ground-based photometric data were then interpreted on the basis of 
simple radiative transfer models . Special attention was given to the case 
of individual, spherical, macroscopic particles for which the scattering 
phase function is proportional to the apparent area illuminated. More 
extreme cases of back-scattering were not considered . Probable ranges 
in the single scattering albedo, and in the general shape of the scattering 
phase function, were defined. Limitations on the mean pex’pendicular 
optical thickness of rings A and B together were also obtained. Results 
indicated that the ring particles are efficient back-scatterers of visual 
radiation. Isotropic scattering is not permitted . Based on an infinite 
optical thickness for the rings, a minimum single scattering albedo 
~0.75 was found. Use of conservative scattering led to a minimum optical 
thickness ~G. 7. 


Using his elementary radiative transfer model, Price (1976) 
recently analyzed the available infrared (10 /im and 20 /im) thermal 
•emission data for Saturn's rings. The ring system was treated as 
an isothermal slab of finite optical thickness , To obtain sufficient 
heating of the ring material, optical single scattering albedos less 
than 0.75 were required. Moreover, the rings could not be isothermal 
To reconcile the results with the analysis of the photometric function, 
the optical back-scattering efficiency of each individual particle 
must be even higher than the maximum adopted . Higher back- 
scattering efficiencies permit lower single scattering albedos, and 
smaller optical thicknesses . 

Further theoretical analysis of the visual photometric function 
data obtained by Price (1973, 1974) is described in this paper. Em- 
phasis is placed on more adequately defining the constraints both on 
the single scattering albedo, and on. the shape of the scattering phase 
function. Extreme degrees of anisotropy are examined. Scattering 
phase functions ranging from total back- scattering to isotropic scat- 
tering are considered . . ■ 
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n. THEORY 

Rrevious studies (Price, 1973.; 1974) have indicated- that primary 
scattering dominates the visual photometric function of Saturn's rings. 
Our analysis will therefore be based principally on the theory of 
“primary scattering, with multiple scattering treated as a perturba- 
tion only. Our radiative transfer model once again assumes the ring 
system to be a plane-parallel layer of particles with an optical thick- 
ness, T^. Each ring particle is described by its single scattering 
albedo, 76, and by its scattering phase fuiiction, p ( (H)), where 
denotes the scattering angle. Once again mutual shadowing is taken 
to he negligible. Interparticle separations are considered to be 
sufficiently large that the angular distribution of radiation scattered 
by each particle is reliably described by the far field approximation. 

Our analysis wiU utilize the entire visual phase curve measured when 
the rings appear wide open with respect to the Sun and Earth. Of the 
available photometric data, these are by far the most reliable. Special 
attention will be given to diffuse reflection near zero phase angle. For 
the geometrical situation of interest, the solar and terrestrial elevation 
angles "above" the ring plane will be taken equal at 26?230 (Price, 1974) 

Writing I, the total specific intensity of radiation diffusely scat- 
tered by the rings , and 7 tF, the solar flux at Saturn, we have (cf: 
Chandrasekhar, 1960) 


^ K’ ® < ’’o' =1 “ ■ 7 • P'® ^ • W 

The parameter denotes the cosine of the angle of incidence, or 
reflection, with respect to the outward normal to the ring plane; its 
ralue will he taken as 0,44124, cf: Price (1974) . The f parameter 
denotes the fraction of the emergent radiation which has suffered only 
a single scattering within the ring system. Equation (1) is the basis 
of our analysis . In principle, since ~ values can be obtained 

directly from photometry of the rings, tight constraints can be placed on 
the permissible combinations of radiative transfer parameters. Before 
doing so, however, two problems must be addressed. First the scat- 
tering phase function must be defined. Second, the primary- to- total 
scattering ratio must be adequately treated. 

Selection of the scattering phase function is based on its meeting ' 
two requirements. First, it should encompass all extremes. Second, 
its general shape should be described by a single parameter . The 
Henyey-Greenstein (1941) function 

( 2 ) 

(1 +g^ - 2g cos (H))V^ 
is the natural choice. General shape of the function is described by the 
g parameter. Extremes include isotropic scattering (g = 01, complete 
forward scattering (g = +1), and complete back-scattering {g = -1).. In 
our analysis, the H-G function will be substituted for p((§)) in equation (1) . 



Since the primary-tp-total scattering ratio is a function of the 
single scattering albedo, of the optical thicfeness, of the scattering 
phase function, and of the geometries of solar illumination and of 
ground-based observation, its evaluation requires careful consideration. 
In general, accurately determining the prianary-to-total scattering ratio 
in a radiative transfer problem which involves an arbitrary Henyey- 
Greenstein phase function can be done only numerically. Straight- 
forward application of the van de Hulst similarity principles is precluded 
both by the situation geometry and by the relative insignificance of 
multiple scattering. Even so, for specific g-values, useful limiting 
f-values can be obtained from known analytical solutions of radiative 
transfer problems which involve similar, but simpler, scattering phase 
functions. From physical considerations, similar phase functions must 
necessarily lead to similar radiative transfer solutions . Our basic 
assumption will be that the corresponding solutions are in fact identical . 
Since multiple scattering is cast in the role of a perturbation only, little 
uncertainty is introduced to the analysis by the approximation. 

Both upper and lower limiting f- mines are used in our analysis . 
Deriving an upper limit is very simple. Obviously, for complete back- 
scattering (g =: -1) only primary scattering can occur, i.e. , f equals 
unity. Deriving a lower limit is not so straightforward, however . From 
Price (1976), the ring particles must back-scatter visual radiation at 
least as efficiently as macroscopic Lambertian spheres. For the latter 
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case, Price (1974) found that the far-field scattering phase function 
could be approximated by 

p(@) = 23 (1 -H X cos {§)) (3) 

•where the anisotropy parameter, x, is set equal to -1. Fig. 1 compares 
scattering phase functions derived from equations (2) and (3) . Isotropic 
scattering is described exactly by setting either g or x equal to zero. 
Moderate back- scattering, represented by x equal -1, is well approxi- 
mated by the Henyey-Greenstein function with g equal to -0.25. Radia- 
tive transfer solutions for x equal -1, given by Chandrasekhar (1960), 
are adopted for g equal -0.25. Using methods employed by Price (1974), 
primary-to-total scattering ratios can be calculated for a variety of 
23, X, T combinations. Results for the geometrical situation of interest 
are given in Table I. All values given are extreme lower f-limits ob- 
tained by adopting the maximum permissible value for the single scat- 
tering albedo from Price (1976). Table I provides information essential 
for the analysis. 


TABLE I 

PRIMARY/TQTAL SCATTERING RATIOS 

(Zero Phase Angle; u = 4 = 0.44124) 

o 


23 

II 

0 

X = -1 


■B 

■1 

1 

7 = C3 

Bfl 


H9B 

B 

.75 

.879 

.690 

.614 

.565 

.955 

.850 

.796 

.759 
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m. ANALYSIS 

Determining the allowable tS-g envelope for the ring particles 

requires a basic knowledge of the visual photometric function for the 

geometrical situation of interest. In rediscussing the observational 

data, Price (1974) showed that the relevant mean value, for rings A 

J? 

I 

and B together, must lie in the range 0.25 s: 0.4, Knowledge of 
the mean optical thickness is important too. For a lower limit, the 
value 0.5 will be adopted. Smaller optical thicknesses would imply 
that stars, and Saturn’s disk, would be readily visible through Ring B, 
a situation contrary to experience. Most probably, the mean optical 
thickness lies in the range 0.5 ^ r ^2. Our present analysis will in- 
clude optical thicknesses in the range 0.5 is r ^ Following Price 
(1976), only IS- values less than 0.75 will he considered. No attempt is 
made to infer optical thickness from our analysis. 

Once are selected, the permitted oi-g envelope can be 

readily calculated using equation (1) , Zero phase angle will be adopted. 
Given g, the only weakness in the precise calculation of oj lies in un- 
certainties in the two parameters, ^ and f , Results of our analysis 
are plotted in Fig. 2. For equal 0.5, maximum ai- values were ob- 
tained by selecting 4^ equal 0.4; for eaual infinity, minimum ui-values 
were obtained by selecting i equal 0.25. In each case, the primary-to- 
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total scattering ratio adopted was relevant both to the optical thickness 
selected and to the scattering phase function characterized by x equal -1, 
c£ ‘jL'able I. For extreme back- scattering (g ~ -1), the adopted f- value 
will be ~ 15-24 percent too small. For -0.25 ^ g ^ -1, the derived 
■w-values will be slightly underestimated; for r equal 0.5, from zero 
to 15 percent; for r equal infinity, from zero to 24 percent. Our final 
conclusions will not be adversely affected, however. Evidently, the 
ring particles must indeed be pronounced back-scatter ers of visual 
radiation, especially if their albedo is low. 
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IV. DISCUSSION 

Current knowledge (Pollack, 1975) suggests that the ring particles 
are macroscopic pieces of H 20 -ice, on the order of a few centimeters 
in diameter. Since knowledge of the surface texture, and general 
shape, of the ring particles does not exist, we will assume for con- 
venience that they are smooth and spherical. Determining if the single 
scattering albedos and scattering phase functions characteristic of 
such particles are permitted by our analysis is important. 

Recently, Veverka (1973) has studied the photometric properties 
of snow and of snow-covered planets. His results are directly applicable 
to the Saturn ring problem. Using measurements of natural snow sur- 
faces at phase angles up to 80 degrees, Veverka concluded that in 
general snow is a quasi-Lambert scatterer, with a characteristic 
Minnaert constant, k, in the range 1.04 1.55. Even in an extreme 

case of specular reflection (a "glazed rain mist"), k was less than 2 
for the phase angles considered. In evaluating the scattering properties 
of snow-covered planets, Veverka derived Bond albedos, Ag, in the 
range 0.51 ^ A^^ s 0.82, with the most likely case A^ >0,67 

For snowballs, analytical expressions for the far field scattering 
phase function were derived for extreme values of the Minnaert con- 
stant, k (viz. , k = 1 and k = 2). Formulae given by Horak (1950) were 
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used. Specifically, we have 

p(@ ) |sin ® . cos^ ® - ® . cos ® + siu^ ©I , k= 1 (4) 

and 

p(d))=|{l-cos®}® , k=2 . (5) 

Both scattering phase functions are shown in Fig. 3. For comparison, 
relevant Henyey-Greenstein functions are also shown. Evidently, the 
scattering phase function for snowballs can be adequately described 
by a Henyey-Greenstein function if the anisotropy parameter lies in 
the range -0.30 s g s: -0.35. Snowballs are superimposed on Fig. 2, 
Low albedos, with very pronounced back-scattering, are also 
permitted by the analysis. Examples of astronomical bodies known 
to have such scattering properties are the Moon and Mercury. Their 
photometric properties are almost identical. Like information for 
other astronomical bodies is not yet complete. Data on the Bond 
albedo of the Moon and on the lunar phase function are readily available 
(Allen, 1963). The lunar albedo is 0.068. The lunar phase function, 
normalized to unity when integrated over all solid angles, is shown 
in Fig. 4. Shadowing within the lunar surface layers is responsible 
for the pronounced back- scattering directivity. For comparison, 
relevant Henyey-Greenstein phase functions are also shown. The 
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lunaT phase function can be adequately described by the Henyey- 
Greenstein phase function if the anisotropy parameter lies in the range 
"0.55 ^ g ^ -0.5, The lunar parameters are also superimposed on Fig. 2. 

Evidently, ring particles TS'ith lunar-like scattering properties 
are excluded by the analysis. Their single scattering albedos would 
be too small for the permitted back-scattering efficiency. Most likely, 
an increase in reflectivity would cause a corresponding decrease in 
back-scattering directivity (Oetking, 1966), Even if the lunar albedo 
were increased by a factor 3 to reproduce the probable values for the 
most reflective asteroids (Chapman, 19Y6), reconciliation with the 
analysis would still be extremely difficult. It might even be impossible. 

More typical, perhaps, of a ring particle might be an average 
grain of lunar surface dust. Their relevsnt photometric parameters 
would be w ~0. 18, g ~-0.25 (Hapke, 1963). Compared with the lunar 
surface as a whole, much less back- scattering directivity occurs 
since surface shadowing is negligible. The lunar dust parameters are 
plotted in Fig. 2. Evidently, such particles are irreconcilable with 
our photometric analysis of the rings. 

Interestingly enough, snowballs encounter no such difficulty. 

Pure spherical particles are in excellent agreement with the analysis. 
Significant impurities in the snow can lead to a discrepancy, 
however. Veverka (1973) has discussed qualitatively their effects. 
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Two kinds of contaminants are possible. First, impurities which do 
not hinder multiple scattering within the surface snow layer to any 
appreciable degree may be present. Examples are frosts of COg, 

CH^, NHg, etc. Second, impurities which effectively prevent multiple 
scattering may be present. These will radically alter the photometric 
properties of a snowball- Most mineral contaminants fall into the 
latter category. Introduction of significant mineral contaminants 
should significantly decrease the single scattering albedo of snow- 
balls, while increasing their back- scattering directivity. Both effects 
will cause movement of the '’snowball” 'S-g envelope in Fig. 2 towards 
smaller ^-values and more negative g-values. 

Quantitative studies of the photometric properties of dirty snow 
have been carried out by Veverka (1976). Contamination was simulated 
by mixing powdered charcoal with magnesium oxide. Geometrical 
albedos, and phase integrals, relevant to smooth planetary bodies, at 
X 5400A, have been evaluated for contamination varying from zero to 
lOG percent by mass. Two data points are superimposed on Fig. 2. 
Evidently, very little mineral contamination of the ring particles is 
permitted by our analysis. Less than ~5 percent contamination would 
seem to be an upper limit . 

Lebofsky, Johnson, and McCord (1970) have studied the visual 
reflection spectrum of Saturn’s rings in detail- Either pure water 
frost, or silicate, particles are inconsistent with their measurements. 
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Part frost-covered silicates and mixtures of frosts with other com- 
pounds, perhaps modified by ultraviolet or high-energy particles, 
were considered to be possible ring constituents. Our photometric 
-.analysis is compatible with their conclusions. 
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Special phase functions, for which analytical radiative 
transfer solutions exist, are described by the anisotropy 
parameter, x. Relevant Henyey-Greenstein functions, with 
anisotropy parameter, g, are shown for comparison. 

Combinations of single scattering albedo, 75, and anisotropy 
parameter, g, permitted by the visual photometric function 
for Saturn’s rings, for optical thicknesses in the range 
0.5 ^ T ^ ®. Regions labelled N/A are excluded by the 
analysis. For comparison, the scattering parameters of 
the Moon, of lunar dust, and of H 2 O snowballs, are also plotted. 
Simulated ’’dirty” snowballs, composed of mixtures of magnesium 
oxide and charcoal, are shown for comparison (Veverka, 1976). 
The MgO/C ratios refer to weight, not volume. 

Scattering phase functions for snowballs, characterized by 
the Minnaert constant, k, equal 1 and 2, are compared 
with relevant Henyey-Greenstein functions, with anisotropy 
parameter, g. 

Lunar phase function compared with relevant Henyey- 
Greenstein functions with anisotropy parameter, g. 
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ABSTRACT 

Harris (1975) has suggested that the maximum size of particles 
in a planetary ring is controlled by collisional fragmentation rather 
than fay tidal stress. While this conclusion is probably true, esti- 
mated radius limits must be revised upwards from Harris* values of 
-a few kilometers by at least an order of magnitude. Accretion of 
■particles within Roche's limit is also possible- These. considerations 
affect theories concerning the evolution of Saturn's rings, of the 
moon, and of possible former satellites of Mercury and Venus. In 
the case of Saturn’s rings, comparison of various theoretical 
scenarios with available observational evidence suggests that the 
rings formed from break-up of larger particles rather than from 
ofTginal condensation as small particles. 
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INTRODUCTION 


Harris (1975) has argued that particles in a ring around a 
planet would rapidly be reduced by collisions to a size much smaller 
than the maximum size limit imposed by tidal stress (Jeffreys, 1947). 
Jeffreys' tide-limited radius was about 100 km both for icy bodies 
in Saturn’s rings and for rocky material near the earth's surface, 
but Harris finds an upper limit due to collisions of -1 km for ice 
and “3 km for rock, apparently independent of the size and distance 
of the planet. Although Harris* work is valuable in calling atten- 
tion to colllslonal effects in rings, there are a number of weak 
points in his discussion. In this paper we present new esti- 
mates for the size limits and discuss uncertainties in the 
cahiulations. In addition we discuss the size distributions 
which might evolve from various initial conditions and the implica- 
tions of these results in theories of planetary evolution. 
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FRACTURE CRITERIA 


Harris noted that a body which receives upon impact an energy 
per unit volume greater than the maximum storable elastic energy 
must fracture. He compared this energy limit to the kinetic energy 
per unit volume provided by a collision, and then solved for a value 
of the collision velocity above which fracture would be certain. 
However, we have identified two necessary corrections to this calcu- 
lation, each of which raises the limiting value of the collision 
velocity: First, the energy per volume prc 'ided by the kinetic energy 

before the collision is at least four times smaller' than the 
expression used by Harris. Second, only a small portion of this 
energy Is available for fracture. We shall consider each of these 
corrections in turn. 

The kinetic energy which goes into a collision of two masses, 
ra and H, is 


1 niM 

2 (m+H) 


0 ) 


where v is the relative velocity before collision. Eq. (l) repre- 
sents the kinetic energy before collision less the minimum kinetic 
energy required after collision to conserve momentum. Another way 
to define this expression is the kinetic energy of the system before 
collision in a center-of-mass -eference frame. Dividing eq. (l) by 
the volume of both particles* we obtain 



B-6 


E/ volume = pv^ 


mH 


Cm+H) 


( 2 ) 


where p is the common density of the particles. However, Harris* 
expression for the kinetic energy per wnlt volume was ipv^, at least 
four times the correct value given by eq. (2). Therefore Harris* 
values for the upper limit to non-fracturing velocities should be at 
least twice the values he presented and perhaps an order of magni- 
tude higher for typical collisions of non-equal masses. 

Why does Harris’ expression differ from ours? He assumed that 
all the strain energy went into the smaller body, rn. However, we, as 
well as Harris, are considering criteria for destruction of the 
larger bodies, and this extreme assumption would make such destruc- 
tion Impossible at any velocity. Even If this assumption is made, 
the kinetic energy density is -jpv^H/ (nrJ-M) , which still differs sub- 
stantially from Harris' expression unless tn « H. 

The second correction involves Harris' assumption that all of 
the kinetic efiergy before collision is available for straining the 
material. Under this assumption fracture occurs if the energy 
exceeds the maximum storable elastic energy, where S Ts 

the yield strength of the material and Y Is the elastic modulus. 

Thus in his model Harris requires fracture when the collision energy 
exceeds 5^10^erg/cm^ for rock and for ice. However, a 

great deal of the collision energy goes into rotational or trans- 
lational kinetic energy of the particles as they burst apart in 
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fracture. Some energy is lost as heat. Moreover, preferential 
dissipation of energy near the contact site might allow only local 
chipping while preventing the widespread stress needed for cata- 
strophic fracture of the larger body. 

■Such factors may explain why experimental collisions have 
required more than the maximum storable elastic energy to yield 
significant fragmentation, as pointed out by Hartmann ^ (1975) 

In a footnote. Moore and Gault (1365) found that, for basalt 
particles with m/M = 8.4x10“^, a coll Ts Ion velocity of km/sec 
can chip only -31 of the mass from the larger body. The energy den- 
sity in that case is 2. 3><l0^erg/cm^. Moore and Gault obtained 
catastrophic collision when they increased the velocity to 2 km/sec 
(S^S-IO^erg/cm^) . " We have performed experiments at lower Velocities 
that corroborate these results. We find that 1,4 gram bullets fired 
at “340 m/sec generally shatter various igneous and coherent conglom- 
erate rock targets only if the target mass is <50 times the bullet 
mass (Fig. l). This result corresponds to a required energy of 
-3"10^erg/cm^, in excellent agreement with Moore and Gault. For 
comparison we note that Gault and Wedekind (I969) obtained -2xl0^erg/cm^ 
as the energy density required for shattering glass targets. These 


"We define "catastrophic collision" or “shattering" as any case in 
which the largest fragment after collision is less than half the 
mass of the largest mass before collision. In irast collision experi- 
ments there Is a sharp transition over a narrow range of energy 
densities from minute local chipping to massive body fracture. 


EEPRODUCIBiLTry OF THE 
PAGE IS POOR 
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results show that substantial fragmentation requires -50 times the 
energy that Harris considered sufficient for fracture of rock (Table l) . 
Therefore, an additional increase of nearly an order of magnitude is 
indicated for estimates of velocities sufficient to yield fracture. 

For collisions of ice particles, it is reasonable to assume 
that the required energy for fragmentation is also -50 times Harris* 
estimate, or 2.5>«l0^erg/cm^. Relevant experimental data on ice are 
rare, Hartmann (publication in preparation) has performed experi- 
ments in which water ice particles impacted a much more massive 
surface. He finds catastrophic fragmentation at velocities above 
-700 cm/sec. If we assume that all the collision energy goes into 
straining the small projectiles we obtain -3’<10^ergs/cm^ for the 
critical energy. Since some unknown amount of the energy goes into 
the large target, the actual energy per volume required would be 
somewhat less than this value. For example, if half the energy goes 
into the target, the critical energy density is -1 ,5>^10^ergs/cm^. 

Harris' value, ergs/cm^, would only be consistent with this 

experiment If -S 8 ^ of the energy goes into the target. We have also 
fired bullets at 3^0 m/sec into hard-packed snow spheres (p - 0 - 5 ). 

We find (Fig. 1) that fragmentation occurs for target masses <10^ 
times the bullet mass, i.e. for energies greater than ' 3 =< 10 ®erg/cm^ , 
in good agreement with our estimates for ice. 

Estimates of the critical energy per volume required for frag- 


mentation are summarized in Table I. We emphasize that all these 
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results must be applied with great caution because of uncertainty 
regarding the environment and materials appropriate to collisions 
in planetary rings. For example, fracture may occur at lower 
collision energies If large bodies contain pre-existing flaws or 
are otherwise weaker than the samples of rock and ice we describe. 
Moreover, we have no mechanical data on material at the low tem- 
peratures which might be appropriate in a ring around an outer 
planet. 
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PARTICLE SIZE LIMIT 

Harris equated his limiting value of fracture velocity with the 
surface ‘'escape'' velocity for a ring particle of radius R: 


• V = 



(3) 


He assumed that typical collisions take place at velocities greater 
than this escape velocity, so eq. (3) permits solution for a value 
of R above which fracture must occur, inserting Harris' own 
values for limiting v and for assumed p, we obtain R = 4 km, for 
rock and 1.3 km for ice. The corrected value of the limiting v, 
discussed above, would make R many times larger. 


One might further question Harris' use of the two-body escape 
velocity in a case where both bodies are inside or near the Roche 
limit of a much more massive third body, the planet. The colliding 
bodies' spheres of influence have sizes comparable to the bodies' 
physical sizes.. The two-body escape velocity Is- irrelevant In this 
case. Though Harris notes that collisions inside the Roche limit 
could occur with arbitrarily small velocities, he nevertheless used 
the two-body escape velocity. 
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We have performed numerical experiments using the restricted 
three-body problem to test Harris’ approximation. The motion of an 
infinitesimal particle under the gravitational influence of a planet 
and of a finite ring particle In circular orbit was followed. The 
primary was given Saturn's mass and the finite ring particle was 
placed at distances corresponding to various points in Saturn's 
rings. The infinitesimal particle was started in various circular 
orbits about the planet and was followed until impact with the 
finite ring particle. The experiments were repeated for various 
masses and radii for the ring particle. (An example Is shown In 
Fig. 2 .) Also experiments were performed v/ith the earth as primary 
and with the finite particle Inside Roche's limit- Collision velo- 
cities were generally somewhat less than, but in each case within 
-25^ of, the "escape" velocity used by Harris. Numerical experiments 
by Reta Beebe (private communication) appear to support this conclusion. 
Thus Harris’ approximation seems justified for an order-of-magni tude 
estimate. Our numerical experiments with the three-body problem 
required one body in the collision to be infinitesimal. We can only 
assume that our collision velocities are appropriate for the case of 
comparable masses. 

Accordingly, we may use eq. ( 3 ) to evaluate the upper limit to 
the size of collision survivors, for a given limiting collision 
velocity, v, as discussed in the previous section. We obtain v 
by equating the energy per volume provided by a collision, eq. (2), 
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with our estimated requirement for fracture (Table I}, 3^10^erg/cm^ 
and 2xl0^erg/cm^ for rock and ice, respectively. If we assume that 
the largest particles have frequent collisions with others of nearly 
the same size, we may sat m = M in eq. (2) to obtain maximum radii 
of '*20 km for ice and ~70 km for rock. These values are 20 times 
Harris' upper limits for the reasons already discussed. The uncer- 
tainties in adopted critical energy for fragmentation do not affect 
this calculation significantly as the velocities, and hence the 
radii, go as the square root of this adopted quantity. 
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PARTICLE SIZE DISTRIBUTION 

The size limit derived in the previous section depended on the 
assumption that the largest particles have sufficiently frequent 
collisions with others of comparable size. This assumption is not 
realistic for two reasons. First, any real collisional history, 
even if it starts with only large bodies, would probably evolve 
toward some approximation of a power law distribution at the expense 
of the large particles (Hartmann 1969). Eventually there would be 
so few particles of the largest size that collisions among them- 
selves would be rare. Second, collisions Tn a planetary ring would 
tend to damp excess relative velocities, yielding nearly circular 
orbits throughout most of the ring. Circular orbits would permit 
th^’large objects to coexist; there would be few mutual collisions. 
(Regions of orbital resonance within the rings enhance orbital 
eccentricities, but these regions are so small as to contain only a 
small portion of the large particles, if any.) Because of these 
considerations the size limit derived In the previous section Is not 
strict; some larger particles might survive. However, we do expect 
collisions to keep the largest particles somewhat smaller than 
Jeffreys* tidal limit because at that limit a particle is, by definition, 
on the verge of being torn apart by the tidal force; virtually any 
collision would be catastrophic. 



ColUsTons might even permit the larger particles to grow 

toward the tidal limit by accretion. A small particle striking a 

larger one might' lose sufficient kinetic energy to prevent escape 

from the gravitational field of the larger one. Accretion might 

seem to be an unlikely process in a planetary ring, such as Saturn's, 

which lies inside Roche's limit where the self-gravitation of a 

particle of non-cohesive material with p = 1 is insufficient to keep 

it from being torn apart by tides. However, consider a small 

particle resting on the surface of a cohesive, much-larger par*" ; 

of radius R smaller than Jeffreys' limit. The tidal acceleration 

tending to pull the small particle off the surface is at most 2 GM^R/r 

where M and r are the mass and distance of the planet. The gravi- 
P 

tational acceleration of the smaller particle toward the larger one 
is GM/R^. For Saturn the binding force dominates for r greater than 
about 1.15 Saturn radii, which is Just inside the inner edge of 
ring C. The surprising result that accretion is possible throughout 
Saturn's ring system is obtained because the maximum tidal strain 
occurs in the larger particle where it can be resisted by cohesion, 
rather than at the interface of the two particles." Whether or not 
accretion is efficient and actually dominates over erosional pro- 
cesses requires a more detailed model of particle interactions, 
Including effects of particle rotation. 


'In our calculation one particle Is much smaller than the other. As 
the size ratio Increases toward unity, tidal acceleration dominates 
at greater distances from the planet. Hence v/e might expect the 
texture of regoUths on ring particles to become finer closer to 
the planet. 
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Given oui* understanding of the nature of collislonal processes 
in rings, what can we say about the particle size distribution to be 
e:<pected from various initial conditions? First we consider 
evolution from a population of large bodies. Suppose a satellite 
or a passing sun■^o^bite^ comes close enough to a planet to be dis- 
rupted by tides. Pollack (1975) essentially ruled out the possibility 
of a large satellite moving into the region of tidal disruption by 
Saturn. He noted that tidal evolution of a satellite's orbit would 
cause it to spiral outward from the planet unless it was already 
inside the synchronous orbit distance v/hich Is well within the ring 
system. However, that statement is only correct for low inclination 
orbits; a highly inclined orbit could contract into the region of 
tidal break-up (Greenberg 197^). A satellite might also enter this 
region after a near-collision with another satellite. 

A large body which enters the region of tidal disruption is not 
shattered as in grinding or collisional break-up into a power law 
spectrum dominated by small particles; rather, it is pulled apart 
into pieces just under Jeffreys' size limit. This is the picture 
described by Harris and by Jeffreys and it seems reasonable to us. 

Even in this first stage of breakup some smaller bodies might be 
created, but pieces near Jeffreys' limit would be expected to dominate 
the population. 
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Elhe expression "Jeffreys^ limit" Is used here to describe the 
maximum size of particles able to resist tidal disruption, it does 
not refer to any particular numerical value, because there is a 
great deal of uncertainty in any evaluation. Jeffreys calculated 
his limit of TOO km using a value of 3^10®dynes/cm^ for the tensile 
strength of ice. This strength was based on studies of ice cliffs 
in nature, Jeffreys recognized that laboratory results gave some- 
what higher tensile strengths. A better estimate for laboratory 
strengths might be -3>=10^dynes/cm^, (cf. Johnson and McGetchin ‘'973). 
Since the radius limit is proportional to the square root of the 
tensile strength In Jeffreys' formulation, ice particles as large 
as 300 km might conceivably survive tidal disruption at the distance 
from Saturn of the present rings. Rocky particles as large as 
3000 km might survive tidal disruption at the same location because 
of their still greater tensile strength.] 

We now consider the evolution of these bodies subsequent to the 
initial tidal break-up. Since the particles’ orbits all originate 
near the same point, the site of the tidal break-up, the orbits must 
overlap and collisions are inevitable. The collisional time-scale 
is difficult to est imate wi thout ad hoc assumptions about the orbits 
and number of particles, so we shall simply make the reasonable 
assumption that collisions are frequent during the age of the planet, 
initial collisions v;ould predominantly involve particles near 
Jeffreys* limit colliding with one another- As we have shov/n, 
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catastrophic collisional fragmentation must occur- Size distributions 
of collisional fragments would be strongly enhanced at a smaller size 
range as in the power law distributions described by Hartmann (1969). 
Eventually, so few particles near Jeffreys' limit would remain that 
collisions with others of nearly the same size would be unlikely. 
Nevertheless, they could still be catastrophically fragmented by 
collisions with much smaller bodies if the relative velocities had 
not been damped out by this time. In this case, especially if 
original relative velocities are extremely high, "super-catastrophic” 
collisions could yield much smaller maximum fragment sizes then we 
have previously calculated. 

Even if the particles are quickly damped into circular equatorial 
orbits, the last remaining particles near Jeffreys' limit would still 
suffer catastrophic breakup in low energy-density collisions because 
they would be strained nearly to the point of fracture by tides anyway. 
However, such break-up would yield pieces of size approximately half 
Jeffreys' limit rather than power-law fragmentation. In this latter 
case the total population's size distribution would be similar to 
the power-law distribut ions Common for hrstorfes of coiTisionar 
comminution, except that there would be an excess of particles at 
about half Jeffreys' limit. 

Further catastrophic fragmentation of the larger bodies in the 
population cannot take place once the orbits are circularized and 
the number of larger bodies is reduced, because small bodies do not 



have enough energy to disrupt larger ones and larger bodies never 
hit one another. Gradual erosion by small bodies cratering larger 
ones might be important. Pollack et al . (1973) have suggested that 
impacts by meteoroids from outside the Saturnian system may have 
dominated the process and created the population of centimeter-sized 
objects observed In the rings. (See next section.) On the other 
hand, substantial accretion of ring particles onto one another would 
also be plausible and they could grov/ in size. However, the growing 
objects would be poorly consolidated and collisions would readily 
break them up again. Thus a steady state betv/een accretion and frag- 
mentation is conceivable. Finally we note that particles less than 
-1 cm in radius could be removed by radiation effects such as 
Poynti ng-Robertson drag (cf. Pol lack et al . , 1973). 

Another possible starting point for a ring might be a collision 
of two bodies near a planet. In this case the Initial size distri- 
bution Is approximately a power-law. If the largest pieces are 
greater in size than Jeffreys' limit, they will be broken by tides, 
yielding a distribution with an excess of larger particles. Subse- 
quent collisional evolution v/ould be indistinguishable from that of 
a particle tidal ly disrupted by a planet as described above. 

So far we have discussed evolution of rings by tidal and 
collisional break-up of larger bodies. Alternatively, suppose rings 
were initially formed of small particles. Particles might grow by 
condensation of material onto some initial nuclei- We would expect 
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particles produced this way to have a small size range. Pollack 
( 1975 ) estimates theoretically an upper size limit of "lO cm for 
such particles in Saturn's Rings. Accretion, as described above, 
would be impossible with all particles of comparable mass. Frag- 
mentation v;ould also be impossible, since high eccentricities and 
inclinations are unlikely in a ring vjhich condensed from a gas in 
this way. The size distribution would probably remain stagnant. 
Substantial orbital evolution due to radiation and Interaction with 
resonance regions v^ould undoubtedly take place, but this evolution 
is not yet well understood. Enhanced eccentricities due to resonances 
could yield some local collisional grinding. 

In summary, we see no reason to expect collisional evolution of 
a ring of large particles to yield particles exclusively less than 
a few kilometers in radius as Harris claimed. The population probably 
would be dominated by small particles as in the power-law distribution 
typical of collisional processes. On the other hand, a ring which 
formed of only small particles would not be expected to evolve to 
Include any larger ones. 
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APPU CAT IONS 
a. Saturn's Rings 

Jeffreys (l9^7) suggested that, since Saturn's rings are only 
a few kilometers thick, they were not formed by disruption of a 
large satellite. Harris, with his result that collisions would 
grind even the largest particle to less than 1 km in radius, found 
Jeffreys' argument untenable. We now find it plausible that such 
evolution might permit survival of bodies many tens of kilometers in 
radius. Indeed, the observational evidence is entirely compatible 
with the existence of large bodies. 

Pollack (1975) has summarized radar, radio, and infrared data 
that demonstrate that most of the cross-section in the rings Is due 
to particles of centimeter scale. Yet these data do not rule out 
the presence of a substantial population of much larger bodies. In 
fact, some recent data seem to require large bodies. For instance, 
quadrant asymmetries in the azimuthal brightness distribution 
(Rietsema , 1976) suggest that at least the order of 10^ of 

the reflected light comes from synchronously rotating particles, 
although a detailed model has yet to be developed. The most plausible 
process that would yield synchronous rotation 1s the raising of tides 
on the ring particles by Saturn, but this process Vi/orks on tlme- 
scaTes shorter than the age of the solar system only for particles 
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>50 m in radius. Moreover, recent observations of local brightness 
variations within the rings by Michael Price and Otto Franz arc 
tentatively interpreted by Price (private communication) as due to 
statistical clumping of particles at least 100 km in radius. 

Although further work is necessary to demonstrate beyond doubt 
the existence o'^ large particles in Saturn's rings, these data seem 
to rule out initial condensation of the rings as small particles. 

Indeed, ve can show that the observations are consistent with a 
power-law size distribution with an ejtponent typical of fragmental 
populations. We may attempt to constrain a power-law distribution 
to the observations that (a) most of tSie cross-section is in 
particles of order 1 cm and (b) -lO^ of the cross-section is in 
particles larger than 50 m radius. We truncate the distribution 
below “0.3 cm radius to account for Pcynt ing-Robertson depletion 
(cf. Pollack 1975). We also truncate above -300 km which is roughly 
Jeffreys' limit and because larger particles v^ould have been observed 
when the rings were edge-on. As shown in Fig. 3, a power-law with 
slope of -3.3 roughly satisfies the constraints, if the slope were 
steeper than that, constraint (b) would be violated^ if it were less 
steep, constraint (a) v/ould be violated. This result is significant 
because populations with exponents of -3-0 to -3-5 are typical of 
fragmentation processes at subsonic velocities (Hartmann, 1369). 
Dohnanyi (1969) similarly calculates a value of -3.5 for a general 
asteroidal fragmentation process. 
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The constraints are probably more closely satisfied by steepening 
the distribution at small sizes to about -3*5 and producing an excess 
at larger sizes (two such curves are shown in Fig. 3 by dashed and 
dotted lines). As we have shov/n, an excess of larger particles Is 
to be expected for several evolutionary scenarios, it is thus 
gratifying that the available observations of the rings seem to 
require a f ragmen tal - 11 ke power- law distribution, probably with an 
excess of large bodies with radii of a few to a few hundred km. 

We thus conclude that the data are completely consistent with origin 
by. tidal and colllslonal break-up of one or a few large parent 
bodies. 

b. Terrestrial planets 

Harris used his results to refute Opik’s (1972) theory that 
the Earth's moon formed out of 200 km particles In a protomoon ring. 
Harris noted that 3 km particles would take an unacceptably long time 
to evolve tidal ly out beyond Roche's limit to where they could accu- 
mulate into a moon. However, we have shown that there is no reason to 
require ring particles to be as small as 3 km. Particles only a 
few times larger could have completed the necessary tidal evolution 
in less than j aeon. On the other hand, Spik's estimate that 
particles would spiral out to Roche’s limit in only 10*^ yr is 


probably much too short, because collisions would undoubtedly reduce 
the mean particle size well below 200 km as we have shown. It must 
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be noted that outv/ard evolution of the large particles would be 
slowed by drag due to the small-size portion of the population. 

This drag would slow the evolution of the larger particles, but 
could provide a mechanism for removing the smaller ones to the site 
of lunar accretion. Given uncertainties in the size-distribution 
of collislonal fragmentation, Opik's theory remains viable, 

Harris further applied his results to Burns^ (1973) and Ward 
and Reid's (1973) suggestion of the possible former existence of a 
satellite of Mercury or Venus. If a satellite spiralled in to the 
Roche limit and broke up, Harris suggests that the pieces would have 
been too small to have been removed by continued tidal evolution. 
How'^v^r, many of the fragments could be much larger than Harris' 
upper limit, so former satellites cannot be ruled out. 
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Tab^e 1: Critical energy per volume required for fragmentation 

Critical Value 

Source Material Method (erg/cm^) 


Harris (1975) 

Typical Rock 

i sVy 

5 X 10^ 


Moore £ Gault (1965) 

Basalt 

Experiment 

(2.3 to 5) 

X 10^ 

This report 

Various Rock 

Experiment 

(2 to 5) X 

lo^ 

Gault & Wedekind (1969) 

Glass 

Experiment 

2 X 10^ 


Adopted value 

Rock 

Tentative estimate 

3 X 1Q7 


1 

Harris (1975) 

J ce 

a s2/y)7 

5 X 1Q3 


Hartmann 

ice 

Experiment 

<3 X 10^ 


This report 

Packed Snow 

Experiment 

(2 to h) X 

10^ 

Adooted value 

Icy Material 

Tentative estimate 

2 X 1Q5 
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CAPTIONS 

Fig. 1: Mass of the largest fragment of a rock or snovf target after 
impact with a small projectile plotted as a function of the 
initial mass ratio. The corresponding kinetic energy per 
unit volume v/hlch goes into the collision is shovjn on the 
Upper scale. 

Fig. 2: Representative trajectories for infinitesimal ring particles 

in initially circular orbits near a 100 km particle which 
moves in a circular, coplanar orbit. The coordinate system 
rotates with the orbital angular velocity of the large rlp'i 
particle. Trajectories of small particles approaching on 
interior circular orbits are similar to those shown in the 
figure, but would approach from the lower right. Impact 
points are labeled with the impact velocities In m/sec. 
Impact velocities are comparable to the escape velocity. 

This relation was found to hold true for various radii, 
densities and distances from the planet. 

Fig. 3 ‘ Possible incremental size-frequency distributions for 

Saturn's rings A and B, assuming particle cross-sections 
account for 50^ of the area of the rings. Hachured lines 
indicate truncations discussed in the text. Observations 
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{mply that most cross-section is in size range (a) and 
is in range (b) ; the equivalent percentages for the 
three distributions are indicated, A single power- law 
with slope - 3.3 on this log-log plot provides a fair fit 
(- 3.3 corresponds to a cumulative mass population Index of 
0.77). The dashed and dotted distributions may show a 
better match to the adopted constraints. The total mass of 
the ring system, due mainly to the largest particles in 
all three distributions, are 2x10”^, 2x10**^, and -10“^ for 
the solid, dashed, and dotted curves, respectively. Only 
the latter starts to encroach on upper limits for the mass 
of the rings suggested by Saturn’s coefficient. 




V 

ring particle 

1 6M / CM^ 

■ 

R = 100 KM 

ESCAPE VELOCITY = 75 M / SEC. 

\ 

TO SATUR^ 

\ 

1 - 10® KM 
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ABSTRACT 

Feasibility of using the rotational Raman spectrum to investigate 
the physical structure of outer planet atmospheres is examined. By 
selecting distinct wavebands spaced throughout the visible spectrum 
for observation, a wide range in physical depth can be probed. On the 
basis of a semi-infinite, homogeneous, pure Hg atmospheric model, 
computations of the strengths of the S{0) and S(l) lines have been made 
for a wide range in physical conditions. For each characteristic wave- 
length, the S(0) ; S(l) line strengths' ratio is independent of the presence 
of aerosol particles. It provides a useful estimate of Boltzmann 
temperature. In contrast, the absolute strengths of the S{0) and S(l) 
lines are extremely sensitive to haze. Extremely small quantities of 
aerosol particles can be detected even at high altitudes . Initial applica- 
tions of the Raman probe technique to Jupiter and Uranus are reported. 
For Jupiter, a weak aerosol haze appears to exist at an Hg column 
density < 17 km-amagat. For Uranus, aerosol particles appear to 
be present at an Hg column density < 104 km*amagat. In addition, 
a significant temperature inversion may be present in the same region 
of the atmosphere . 
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1. INTRODUCTION 

In principle, the Raman effect provides a powerful astrophysical 
tool for determining the physical structure of planetary atmospheres in 
general. But, in practice, the Rayleigh and Raman scattering cross- 
sections of the predominant gas must be of essentially the same order 
of magnitude for the Raman effect to be readily detectable. Moreover, 
the atmosphere must be optically deep with aerosol particle scattering 
being relatively insignificant. Fortuitously, is one of the only gases 
for which the former condition is met. Development of the Raman probe 
technique is, therefore, directly relevant to the four outer planets. Both 
Uranus and Neptune are prime candidates since the consensus of present 
evidence indicates that their atmospheres are deep, essentially clear 
above the clouds, and composed of nearly pure H 2 . For wavelengths 
less than about 6000 1, their optical scattering properties appear to be 
determined almost completely by Rayleigh and Raman scattering. 

The Raman effect has two important features which together make 
it effective in probing H 2 planetary atmospheres. First, the efficiency 
of both Rayleigh and Raman scattering increases essentially as the 
inverse fourth power of the incident wavelength. Because the optical 
depth (at X ^ 6000 1) is directly proportional to the effective Rayleigh/ 
Raman scattering cross-section, the mean physical depth to which 
photons penetrate before reflection (i.e., the integrated number of H 2 
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molecules in the line of sight) varies essentially as the fourth power of 
the wavelength. So, the Raman contribution function will move progres- 
sively deeper into the atmosphere with increasing wavelength. Second, 
the intensities of the Stokes and anti-Stokes lines, in the optically thin 
case, are proportional to the populations of their corresponding initial 
states; the intensity ratio is proportional to the population ratio in terms 
of the Boltzmann distribution. More generally, the intensity ratio is a 
function of the population ratio which is related to temperature. 

Utilizing both features, studies of the Raman spectra of both 
Uranus and Neptune over the optical range in wavelength should permit 
coarse definition of the physical structure and temperature of each 
lower atmosphere. Corresponding studies of the ultraviolet spectra of 
both Jupiter and Saturn should provide similar information on their 
upper atmospheres. This paper explores the feasibility of using the 
Hg rotational Raman spectrum to probe the outer planets. 

2. TEE SCATTERING PROBLEM 

Our radiative transfer model consists of a semi-infinite, hazy, 
homogeneous, isothermal, atmosphere in hydrostatic equilibrium. 
Temperatures in the range 0 < T < 150°K, relevant to outer planet atmo- 
spheres, are adopted. It is assumed that the H 2 molecules are completely 
thermalized so that their energy levels are populated according to the 
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Boltzmann distribution. Table I shows the relative populations as a 
function of temperature and energy level. To an excellent approxima- 
tion, all molecules can be taken to be initially in their lowest 
vibrational state (v = 0). For the temperatures of interest, the popu- 
lation of the rotational state J = 2 is negligible. Consequently, all Hg 
molecules can be taken to be initially in the J = 0 and J = 1 states. 

Table n lists the relevant Raman transitions. Wavenumber shifts are 
taken from the adopted tabulation of scattering cross-sections given by 
Ford and Browne {1973). 

Wallace (1972) has derived the equations of radiative transfer 
for the Raman scattering problem. Predictions of the strengths of the 
S(0) and S(l) rotational Raman lines can be reliably based on his approxi- 
mate solution. Using the Wallace (1972) notation, we have 

^ ( 0 ,^ ( 1 ) 

ft i 

where 1^(0, p.) is the specific intensity of radiation emerging from the 
atmosphere which has suffered only Rayleigh scattering, and 12 ( 0 , p) 
is the corresponding intensity of radiation which has suffered only a 
single rotational Raman shift. The parameter p is the cosine of the 
angle of emergence with respect to the outward normal. Physically, 
the parameter is the probability that if a photon colli Jes with a gas 
molecule, it is scattered with no change in energy. The parameter 
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ia the corresponding probability that a photon will experience a change 
in energy relevant to a selected rotational Raman transition. More 
specific definitions of the relevant parameters are 


- FTTkTTX 


( 2 ) 


Note that the definitions of lj^(0,p) and can be expanded to include 
conservative aerosol particle scattering. 

The parameter is the effective Rayleigh-aerosol particle 
scattering cross-section. To avoid adding further to the list of parame- 
ters, aerosol particle scattering is taken to be conservative. Specifically, 


O- = CT 
1 g 




( 4 ) 


where c is the effective Rayleigh scattering cross-section. The 
& 

parameters are the individual mean free paths for Rayleigh 

scattering and for aerosol particle scattering, respectively . Specifically, 





( 5 ) 
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•where N is the number density of Ho molecules, and 
g ^ 

\ ( 6 ) 

"Where N is the number density of aerosol particles, and a is their 
c ^ 

effective cross-section. For homogeneous atmospheres, A. is 

t 

constant with optical depth . 

The parameter is the cross-section for Raman scattering from 
■wavelength 1 to wavelength 2, weighted according to the fraction of H 2 
molecules in the initial state; it is taken as relevant either to the S(0) 
or S(l) transitions. The parameter is the effective cross-section 
for Raman scattering from wavelength 1 to wavelengths other than wave- 
iength 2. It includes all suitably weighted Raman transitions shown in 
Table II, except that covered by X^. Although is dependent on the 
aerosol content of the atmosphere, both X^ and k^^ are not. For a clear 
atmosphere, Fig. 1 illustrates the variation of and with wavelength 
throughout the range 1500 A - 5500 1. Two extreme situations are con- 
sidered. In one, all H 2 molecules are initially in the J = 0 state; in 
the other, they are all in the J = 1 state. 

Determining the contribution which each element of optical depth 
makes to the strength of each rotational Raman line is an essential part* 
of our feasibility study. Any photon escaping from the atmosphere 
which contributes to the strength of either the S(0) or 3(1) lines can 
have experienced only a single Raman shift from its original wavenumber. 
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Its prior, or subsequent, history can have included only Rayleigh 
and/or aerosol particle scattering with no change in wavenumber. Since 
the rotational Raman shifts are small, the basic optical scattering 
parameters remain essentially unchanged . Both the single scattering 
albedo, and the effective scattering mean free path, for Raman shitted 
photons remain essentially identical to those for Rayleigh/aerosol scat- 
tered photons . 

To include the effect of wavenumber shift on optical depth in his 
theoretical treatment of Raman scattering, Wallace (1972) introduced the 
scaling factor, defined by 




or. 4- k. + X. 

Ill 

-f k^ + 


( 7 ) 


The parameter K 2 relates optical depths at wavelength 1 to those at 
wavelength 2 . From the natural vp.riation of the Rayleigh and Raman 
scattering cross-sections with wavelength, all Raman shifts towards 
the red region of the spectrum cause Kg to become less than unity. 
Estimates of Kg can be obtained on the assumption that the Rayleigh 
and Raman scattering cross-sections vary as the fourth power of wave- 
number. If wavenumbers before, and after, a Raman 

shift, we can write 


Kg - 


( 8 ) 
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For the S(0) line, K 2 lies in the range 0.98 ^ ^ 0.92 for 1500l ^ X ^ 

5500 I. For the S{1) line, IC^ can take values of 0.97 ^ K 2 ^ 0.88 for 
1500l ^ 55001. By comparison, Wallace (1972) obtained ~0.95 

for the S(0) line at X 4000 A. Bear in mind that equation (8) refers to 
a clear gas. Introduction of aerosol particles into the atmosphere will 
cause Kg to move back towards unity. Evidently, if the Rayleigh/ aerosol 
scattering contribution function is known, it can also be used for rota- 
tional Raman scattering. 

Analytical solution of the basic transfer equation for a semi- 
infinite, homogeneous, isotropically scattering, planetary atmosphere 
has been obtained by McElroy (1971) on the basis of the two-stream 
approximation to describe the radiation field. His solution has been 
adopted in our feasibility study of the Raman probe technique. If 
lj^(0,ju) is the specific intensity of Rayleigh/aerosol scattered radiation 
diffusely reflected in the /x-direction, we can write 


Il(0,M) 

F 


J 

0 


c(t, 

F 


dr 


(9) 


where ttF is the solar flux at the planet, and c{r,}j.) is the contribution 
function relevant to an optical depth r in the h-direction. The contribution 
function is defined by 


c(t,m) = ^ • J(r) exp (-r/p) 


( 10 ) 
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where J(r) is the source function, given by McElroy (1971) as a 
complicated function of [i, and Mq- In the usual notation, iSq is the 
cosine of the angle of incidence with ref- i' ct to the outward normal. 

For convenience, w- will assume are equal, a valid approximaLion 
for the outer planets . 

Contribution functions for 0.9 ^ oj ^ 1 and 0. 1 ^ /i ^ 1 are shown in 
Fig. 2. Evidently, the major contribution to radiation escapuig from 
the atmosphere comes from optical depths r ^ 2. This is especially 
true if p is small. For a clear atmosphere, Fig. 3 illustrates the vari- 
ation of optical thickness with wavelength for a constant column 
density. Only Rayleigh scattering from H 2 molecules is considered. 

For convenience, all H 2 molecules are taken to be in the v = 0, J = 0 
state. Essentially the same curve would result if all H 2 molecules were 
in the v = 0, J = 1 state. All optical depths, r. have been normalized to 
unity at 5500 1. Plainly, the physical depth probed by the Raman effect 
is strongly dependent on wavelength. Bear in mind, however, that the 
introduction of aerosol parcicles into the atmosphere will progressively 
reduce the wavelength dependence of the physical depth probed . 

Our exploratory investigation of the Raman probe technique will 
consider radiation from an entire planetary disk only. Attempting to 
investigate the variatioi-. of the Raman components across the disk of a 
planet is premature, especially since Uranus and Neptune cannot be 
adequately resolved by ground-based techniques. Predictions of the 
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spectral intensities of the S(0) and S(l) lines will involve geometrical 
ahbedos exclusively. Geometrical albedos, can be readily derived. 

J’or Mq equal, "we can write 

1 1.(0, M) 

p = 2 J — — ii dfi . (11) 

Geometrical albedos, P]^?P 2 refer to the Rayleigh/aerosol particle 
scattering and to the rotational Raman scattering, respectively. 

To calculate P^,P 2 > expressions for 1^(0, jLi)/F and a 1^(0, 

can readily be derived from McElroy (1971). Two-point Gaussian 

quadrature is sufficient to evaluate the integral to 0.1 percent accuracy. 

Determining the accuracies within which the two-stream approxi- 
mation for the radiation field defines the geometrical albedos, Pj>P 2 > 
is important. Accordingly, values were calculated for 0.9 ^ ^ 1.00. 

In Table III, they are compared directly with exact values given by Harris 
(1961). Through the use of the two-stream approximation, pj^ values 
are consistently underestimated by from 11.5 to 14.8 percent. Accuracy 
of the p 2 values can be estimated by applying the finite difference tech- 
nique to the data given in Table III. Results of the exercise are presented 
in Table IV. Evidently, is systematically underestimated by 

17-19 percent. By comparing the strength of the Raman component of 
interest directly with the Rayleigh/aerosol scattering component, 
systematic errors in Pj^.P 2 can be partially canceled. :I'ven so, the 
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Pg/Pl ratio will be underestimated by ~2 - 5 percent. Such errors are 
tolerable in our current feasibility study. 

For convenience, spectral detectability of either the S(0) or S(l) 
line is defined here in terms of the geometrical albedos' ratio, 

In principle, the total geometrical albedo, p,^,, should be used in place 
of pj. For each wavelength of interest, additional photons will have 
been shifted from shorter wavelengths either by single or multiple 
vibrational and/or rotational Raman scattering. Spectral detectability 
(P 2 /P 1 ) therefore be overestimated by the factor p.ji/p^. Bearing in 
mind that the sum of the geometrical albedos resulting from all wave- 
length paths cannot be greater than the geometrical albedo resulting from 
conservative scattering alone, we can use Fig. 1 and Table in to estimate 
an upper limit to P-^/p^j^. From Fig. 1, we have oj > 0.9 for all wavelengths 
in the range 1500 1 - 5500 1. From Table III, p(w = l)/p(2i = 0.9) is 
0.690/0,310. Consequently, spectral detectability cannot be overesti- 
mated by more than a factor 2.2. Since the solar photon flux decreases 
nearly exponentially with wavenumber in the spectral region 1500 1 - 5500 1, 
the factor is most probably much less than two. Indeed, results obtained 
by Wallace (1972) suggest values for P-p/Pj ~1.25 only. 

Theoretical computations were made for both clear and hazy H 2 
atmospheres. Raman contributions to the geometrical albedo are denoted 
p(J = 0) and p(J = 1) for the S(0) and S(l) lines, respectively. For a clear 
atmosphere, the geometrical albedo is denoted 
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aerosol haze, ifc is denoted p. In an obvious notation, the 

relative population (f) of the two initial states is given by 
N{J = 0)/ |N(J = 0) + N{J = 1) f - Results for a clear atmosphere are 
illustrated in Fig. 4. For both the S(0) and S(l) lines, spectral 
detectability (p 2 /pj^) is plotted as a function both of wavelength and 
relative population of the two initial states . For each Raman line, 

, spectral detectability is essentially independent of wavelength in the 
range 1500l - 5500a. It is, however, nearly linearly dependent on the 
population of the relevant initial state. Extremely cold H 2 atmospheres 
will have S(0) Raman components with maximum spectral detectabilities 
~20 percent. For high temperatures (X 150-200'^K), the maximum 
spectral detectability of the S{1) line will be ~13 percent. 

Fig. 4 strongly suggests that the ratio of the S(0) and S{1) 
absolute line strengths might be a useful indicator of temperature for 
those regions of the atmosphere probed by radiation of a selected 
wavelength. But it is essential for p(J = 0)/p(J = 1) to be insensitive to 
the presence of aerosol particles in the atmosphere. To investigate the 
matter, the necessary Raman strength calculations were made for an 
arbitrary wavelength (4000 a), and for an arbitrary situation in which 
equal numbers of H 2 molecules are in the J = 0 and J = i initial states 
(f = 0.5). Calculations were performed for those cases in which the 
scattering mean free path for H 2 molecules {X ) relative to the corre- 
sponding distance for aerosol particles alone (X ) varied in the range 
0 ^ ^ Results are plotted in Fig. 5. Evidently, p(J = 0)/p(J = 1) 
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is insensitive to So, the S(0)/S(1) ratio can indeed be a useful 

to c 

indicator of temperature. 

By contrast, the individual spectral detectabilities of the S(0) 
and S(l) lines are extraordinarily sensitive to the aerosol content of 
the atmosphere. Illustrative calculations of p(J = ®)/P^erosol’ 


p(J = l)/p 


Aerosol 


, were made for an arbitrary wavelength, 4000 A. 


Extreme variations in aerosol scattering were treated by considering 
mean free paths in the range 0 ^ X ^100. Results are plotted in 
Fig. 6. Evidently, the spectral detectabilities of the S(0) and S(l) lines 
are sensitive to the presence of even trace amounts of aerosol particles, 
z.e., for X^X^<2. 


3. APPLICATIONS 
3, 1 Jupiter 

Fast, Poeckert, and Aumann (1974) have detected pure rotational 
Raman scattering from H 2 in Jupiter in the 3500.4 - 3900.4 spectral 
region by means of a photographic correlation technique. The combined 
intensity of the S(0) and S(l} pure rotational lines was found to be 
3.5 + 0.9 (p.e.) percent of the total intensity in this spectral region. 

No information was given concerning which regions of Jupiter were 
studied, so it will be presumed that the Jovian spectra corresponded 
to the mtcgrated light from the entire planetary disk. 



I 


C-16 

To derive their estimate of the combined intensity of the rota- 
tional Hainan lines, Fast et al. selected a priori the relative S(0):S(1) 
intensity ratio. On the assumption that each photon is scattered once, 
the intensities of the pure rotational Raman lines S(0) and S(l) were 
found to be 0.93 and 1.42 percent, respectively, of the Rayleigh intensity. 
The necessary scattering cross-sections were derived from formulae 
given by Dalgarno and Williams (1962) and the polarizability data of 
Victor and Dalgarno (1969) and Victor, Browne, and Dalgarno (1967). 

The required Raman frequency shifts and molecular constants were 
obtained from Stoicheff (1957) . The temperature was arbitrarily taken 
to be 135°K. The S{0):S(1) intensity ratio 0.93:1.42 was used directly 
in the correlation analysis. Obviously, no informacion regarding 
temperature can be obtained from the results. Further, it is by no 
means clear whether the ultimate result of the correlation analysis is 
adversely affected by a priori assumptions regarding the S(0);S(1) rela- 
tive strengths. 

Even so, upper limits to the aerosol content of the relevant 

regions of the Jovian atmosphere can be obtained. For the case m 

point, it can be assumed that all H 2 molecules are in their lowest 

vibrational state. For convenience, ^ H 2 molecules will be placed 

in either the J = 0 or J = 1 state. Fig. 6 can be used to provide an 

adequate estimate of X A for each case. For J = 0, X A = 29"', j: ; 

g’ c ’ g c -12’ 



C-17 


for J = Ij . Evidently, we can write 8 ^ X ^ 

both J = 0 and J ~ 1, the Rayleigh scattering cross-section at 3700 A is 
-27 2 

5.0 X iO cm, . From equation (4), the effective Rayleigh/aerosol 

- 9fi 9 

scattering cross-section lies in the range 4.5 x 10” cm ^ ^ 

-25 2 

2,55 X 10 cm . For unit effective optical depth, the corresponding 
Hg column densicy, must lie in the range 1.5 km-amagat ^ ^ 

8.3 km-amagat. Since the rotational Raman lines are mostly formed at 
effective optical depths less than ~2, significant aerosol particle scat- 
tering must be occurring in the Jovian atmosphere at H 2 depths equiva- 
lent to less than ~ 17 km-amagat. Evidence for aerosol particles so high 
in the Jovian atmosphere was reviewed by Axel (1972). Unpublished 
photometric studies of the Pioneer 10 and 11 imagery by Tomasko (1976) 
confirm the presence of aerosol particles in the regions probed by the 
rotational Raman effect , 

3.2 Uranus 

Belton, Wallace, and Price (1973) detected the Raman effect in 
the spectrum of Uranus between 3880 1 and 4140 1. The spectrum 
corresponded to the integrated light from the planetary disk. Stokes' 
shifts corresponding to both the S(0) and S(l) rotation-i', transitions in 
H 2 are present. A composite Uranus/Moon ratio spectrum comprised 
the observational material. To compare theoretical predictions with 
observation, a synthetic Uranus/Sun ratio spectrum was constructed. 
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Computations by Wallace (1972) of the relative strengths of the Rayleigh- 
T?g man-.qhi fieri components of the spectrum of a semi- infinite, pure 
H 2 , planetary atmosphere were used as the basis for modelling the 
Uranus atmosphere. Both S(0) and S(l) Raman components were included 
with the restriction that the sum of their strengths should equal the 
Wallace (1972) result for the S{0) transition alone. Contributions from 
the higher rotational levels were expected to be small and were neglected, 
Raman components arising from the vibrational transition v = 0 to v = 1 
and from multiple rotational and/or vibrational shifts (cf: Wallace, 1972) 
were included based on the approximation that the shape of the solar 
spectrum below 3850.1 could be represented by a black body at 5800°K. 

The relative S(0);S(1) Raman contributions were 0,3:0. 7. The uncer- 
tainty in their ratio appeared to be between !t0. 1 and iO.2. The fractional 
contribution from the 3(0) Raman transition to the total spectral intensity 
was ~0.03; for the 3(1) transition, it was -^0.08. Although the results 
were consistent with the Belton, McElroy, and Price (1971) dear- 
atmosphere model, neither of two alternatives could be ruled out. 

Specifically, Danielson and Wannier (1973) derived a preferred model 
for the atmosphere consisting of a dear Layer of Hg overlying a dense NH^ 
crystal cloud; the H 2 column density was ~500 km-araagat. Belton, Wallace 
and Price (1973) were able to place limitations on the H 2 column density by 
investigating the effect on the Raman intensities of including an opaque 
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Lambertian cloud layer "within the dear-atmosphere model. For unit re- 
flectivity of the cloud deck, the PJaman scattering obser" 7 ations were found 
to be consistent with any Hg column density greater than ~200 km-amagat. 

By comparison, Prinn and Lewis (1973) have argued on thermo- 
djmamic grounds that a significant methane haze should be present high 
in the atmosphere (H 2 column density ^ 119 km-amagat), with a deep 
dense ammonia cloud layer far below (H 2 column densities 370 km-amagat). 
Belton (1973) included isotropic aerosol particle scattering in the 

semi- infinite, homogeneous, clear atmospheric model to investigate the 
idea. From Prinn and Lewis (1973), the atmosphere can be approxi- 
mated by one in which the ratio of the H 2 volume scattering cross-section 
to the aerosol particle volume scattering cross-section could be 0.5 
at 4000 1 if isotropic scattering were operative. With all the H 2 mole- 
cules in the v = 0, J = 1 state, Belton d. (1973) found the ratio of 
the contribution to the spectral intensity from the S(l) Raman component 
to the sum of the Rayleigh and all Raman components to be 0.07 instead 
of 0.10 for no aerosol scattering. They considered this to be not in 
conflict with the observed ratio ~ 0.08 iO.Ol. 

Upper limits to the aerosol content of the relevant regions of 
the Uranus atmosphere can be obtained on the basis of our approximate 
theory. Once again, it can be assumed that all H 2 molecules are in 
their lowest vibrational state. Also, all H 2 molecules will again be 

Ob' lii" 

1 POOH 


EEPBODUCmOT 
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placed either in the J = 0 or in the J = 1 rotational state. Fig. 6 can 
be used to provide an adequate estimate of X for each case . For 

j = 0, Xg/x^ = ~ ^c/^c ~ ^'^-0*5 * 

Tvrite X ^/X^ > 1. For both J = 0 and J = 1, the Rayleigh scattering cross- 

"2T 2 

section at 4000 1 is 3.6 x 10 cm . From equation (4), the effective 

-27 

Rayleigh/aerosol scattering cross-section is greater than 7.2 x 10 
2 

■ cm . For unit effective optical depth, the corresponding column 

density, N,j,, must be less than 52 km-amagat. Since the rotational Raman 
lines are mostly formed at effective optical depths less than~2, the 
regions of the Uranus atmosphere probed correspond to H 2 column densi- 
ties less than ~104 km-amagat. Within this region, the aerosol volume 
scattering coefficient is at least comparable with the H 2 volume scattering 
coefficient, in agreement with the predictions of Prinn and Lewis (1973). 

Belton et al. (1973) also attempted to investigate the validity of 
the Trafton (1967) model of the Uranus atmosphere. His pressure- 
temperature relation defines a clear inhomogeneous atmosphere in 
which isotropic scattering instead of Rayleigh was assumed. Computa- 
tions of the corresponding relative strengths of the S(0) and Sfl) Raman 
features were .aade. The calculated S(0)/S(1) ratio 0.67 could be 
compared with the observed ratio 0.4. Bearing in mind the uncertainty 
of ±0.01 to ±0.02 in the observed ratio, theory and observation were 


not considered to be inconsistent. 
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Reasons for the presence in the Uranus atmosphere of a thermal 

0 -2 

inversion of at least 140 K at a pressure ~3 dynes cm were recently 
discussed by Wallace (1975). Fig. 4 and Table 1 may be used to 
estimate a Boltzmann temperature corresponding to the 4000l observa- 
tions of the strengths of the S(0) and S(l) Raman lines. For p(S = 0)/ 

PRayleigh^°-°^’ 0.08, f 0.3 , So, 

0.15 <f <0.3. Interestingly enough, the required temperature is greater 
than 150°K. 


4. DISCUSSION 

Feasibility of using the Raman probe technique to investigate 
the physical structure of H 2 planetary atmospheres has been demon- 
strated. But detailed investigations of the outer planets will require 
progress in the areas of theoretical prediction and of observational 
analysis. Theoretical predictions should take account of multiple 
Raman shifts to better determine the continuum level at each selected 
wavelength. Although current neglect of multiple Raman scattering 
does not adversely affect estimates of temperature, it does cause a 
slight overestimation of the aerosol content of the atmosphere. Effects 
of aerosol particle absorption should also be investigated. Uniqueness of 
the radiative transfer solutions should be investigated using inhomogeneous 
atmospheric models. Bear in mind that the Raman contribution function 



■will be significantly affected by the detailed physical structure of the 
atmosphere . Observational analyses should not require a priori 
assumptions concerning the relative strengths of the S(0) and C'D lines. 
Expeditious accurate measurement of the liam'iin component strengths 
fay correlation with the solar spectrum properly requires maximum 
spectral resolution and photometric accuracy. Ideally, the Raman 
"ghost" images of the Fraunhofer spectrum should not require use of 
correlation techniques for their detection. With further development, 
the Raman probe technique could prove to be a valuable tool for in- 
vestigating Jupiter , Saturn, Uranus, and Neptune. 
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TABLE n 
RAMAN SHIFTS 


Transition 

Wavenumber {cm 

-7 = 0 , J = 0"^v = 0, J = 2 

354.39 

v-0, J=l-»v = 0, J=3 

587.07 

v = 0, J=0-^^v = l, J = 0 

4162.06 

v = 0, J = 0-»v~l, J = 2 

4498.75 

v = 0, J=l-^v = l, J=1 

4156.15 

v = 0, J=l‘+v = i, J=3 

4713.83 
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TABLE m 

GEOMETRICAL ALBEDOS: EXACT & TWO-STREAM 


rw 

(a 

p (Exact) 

p (Two-Stream) 

p (Exact)/p (Two-Stream) 

1.000 

.690 

.601 

1.148 

.999 

.636 

.555 

1.146 

■ .995 

.575 

.503 

1.143 

.990 

.534 

.468 

1.141 

.975 

.462 

.407 

1.135 

.950 

.392 

.348 

1.126 

.900 

.310 

.278 

1.115 










TABLE IV 


FINITE DIFFERENCES FOR EXACT & TWO-STREAM CALCULATIONS 



1-t 

(Exact) 

Apj (Two-Stream) 

Apj(Exact)/Ait)j^ 

Ap^ (Two-Stream)/Aoi j 

Ap ^ (Exa ct)/Ap 1 (Two - Stream) 

1.000 

0.001 

0.054 

0.046 

54.00 

46.00 

1.17 

0.999 

0.004 

0.061 

0.052 

15.25 

13.00 

1.17 

0.995 

0.005 

0.041 

0.035 

8.20 

7.00 ' 

1.17 

0.990 

0.015 

0.072 

0.061 

4.80 

4.07 

1.18 

0.975 

0.025 

0.070 

0.059 

2.80 

2.36 

1.19 

0.950 

0.050 

0.082 

0.070 

■ 

1.64 

1.40 

1.17 

0.900 
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Fig. 1 


Fig. 2 


Fig. 3 


FIGXJKE CAPTIONS 

The Rayleigh and Raman scattering parameters, and 
versus wavelength. The atmosphere is assumed to be 
semi-infinite, clear, isothermal, and composed of pure 
Hg* Two extreme situations are considered. Besides 
being in their lowest vibrational state, a^H 2 molecules are 
initially either in the J = 0 state or in the J = 1 state . 

Specific intensity contribution functions, c(r,/i), versus 
optical depth, r. Four choices of single scattering albedo, 
are considere?^ together with three choices for the 
incident and emergent geometry = jUq). The quantity ttF 
is the solar flux at the planet of interest . Computations refer 
to a semi-infinite, homogeneous, isotropically scattering 
atmosphere. The two- stream approximation is used to 
describe the radiation field. 

Optical depth scaling factor versus wavelength. Only H 2 
Rayleigh scattering is considered. All H 2 molecules are 
assumed to be in thetr lowest vibrational and rotational 
states. All values have been scaled to unit optical depth 
at 5500 i. 
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Fig. 4 Strengths of the S(0) and S(l) rotational Raman lines 

as functions of wavelength and population of the initial state , 

The atmosphere is assumed to be semi- infinite, homogeneous, 
isotropically scattering, and composed of pure Hg. Besides 
being in their lowest vibrational state, all Hg molecules are 
assumed to be initially either in the J = 0 or in the J = 1 state . 

The parameter f indicates the fraction of the molecules 
in the J = 0 state; (1 - f) is the corresponding fraction in the 
J = 1 state. The Raman scattering contribution to the geometrical 
albedo, p, is scaled according to the Rayleigh scattering contribu- 
tion. 


Fig. 5 Relative strengths of the S(0) and S(l) Hg rotational Raman 
lines as a function of aerosol content of the planetary atmo- 
sphere . Line strengths are given in terms of the Raman 
contributions to the geometrical albedo. The atmosphere is 
taken to be semi- infinite, homogeneous, and isotropically 
scattering. The aerosol concentration is given in terms of 

X /A. where X is the Rayleigh scattering mean free path, 
jo ^ S 

and X is the aerosol particle scattering mean free path. Com- 
putations were made for a wavelength 4000 1. The molecules 


were considered to be equally distributed between the J = 0 
and J = 1 initial states. Notice the scale change in the abscissa. 
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Fig. 6 Geometrical albedos for the S(0) and S(l) H 2 rotational 

Haman transitions, relative to those for Rayleigh/aerosol 
scattering, plotted as a function of aerosol content of the 
atmosphere. The atmosphere is taken to be semi-infinite, 
homogeneous, and isotropically scattering. The aerosol 
concentration is given in terms of where is the 

Ho Rayleigh scattering mean free path, and X is the 
aerosol particle scattering mean free path. Computations 
were made for a wavelength 4000 A. Two extreme situations 
are indicated. All Hg molecules are initially either in the 
J = 0 state or in the J = 1 state . 
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ABSTRACT 


Multi-color (X5500A - X7600l), narrow-band (lOOA), area-scanning 
photometry has been used to study the wavelength variation in the optical 
appearance of Uranus. Limb -brightening has been detected in two CH^ 
bands i.e. X6X90l and X7300l. . 
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1. WTRODUCTION 

Belton and Vesculus (1975) have recently reviewed the UV-IR 
appearance of Uranus in the context of current knowledge regarding its 
atmospheric structure. Within a predominantly H 2 -* CH^ atmosphere, 
the tops of a thick NHg cloud layer should exist near the 3-4 bar level; 
there is also a very uncertain possibility of a thin, broken, CH^ cloud 
layer near 300 mbars. But whether in fact the Uranus atmosphere is 
so structured remains very much an open question. Two major questions 
are of paramount interest; 

1. Does the deep-dense cloud layer really exist? Is it in 
fact observable? If so, at what pressure level is the top 
located ? 

2. Is the H 2 - CH^ atmosphere above the NH 2 cloud layer per- 
fectly clear? If not, what is the precise nature of the aerosol 
haze? 

Answers may be obtained by imaging the planet from a fly-by spacecraft. 

But, until a Uranus mission is successfully accomplished, further knowledge 
of the planet must come from Earth- or Satellite-based, observations. 

Belton and Price (1973) have developed a straightforward test of the 
hypothesis that clouds are absent from the visible atmosphere of Uranus . 
Based on a clear, semi-infinite, pure H 2 - CH^ gaseous model, it consists 
of measuring the wavelengths at which the character of the center-to- limb- 
intensity profile changes from limb -darkening to limb-brightening, and 
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vice versa. Limb-brightening occurs in certain wavelength regions 
by virtue of the vertical inhomogeneity in the atmosphere which results 
from pressure-related opacity sources. Where the dominant opacity is 
due to pres sure- induced and pressure-broadened transitions in Hg and 
CH^, the intensity profile of the Uranus disk changes character at a geo- 
metrical albedo of 0.23. Seven changeover wavelengths (+20^ were 
predicted viz., 6120l, 6240l; 6620i, 6710l; 6970i, 7420l; 7600l. All of 
the deep CH^ bands beyond 7600l should be limb-brightened, as should the 
central regions of the CH^ bands at 61901 and 7300l. Fig. 1 summarizes 
the Belton and Price (1973) predictions. Available observations, discussed 
below, were found to be in qualitative agreement with the model. But, they 
were not yet numerous enough, nor sufficiently accurate, nor made in the 
critical spectral regions to be conclusive. 

For visual wavelengths (3800l - 5800l), Danielson, Tomasko, and 
Savage (1972) obtained a high-resolution composite image of Uranus. By 
using Stratoscope H, they were able to achieve a Gaussian point spread 

ft 

function with a halt-intensity width 0.2 arc. Except for the presence of 
pronounced symmetrical limb-darkening, no certain disk, markings were 
visible. Any faint belts parallel to the rotational equator must have a 
maximum contrast of 5 percent. The measured limb-darkening did not 
agree with either a deep Rayleigh atmosphere or with clouds high in the 
atmosphere; a cloud deck underlying a finite Rayleigh atmosphere seemed 
to be indicated. Light and Danielson (1973) found later that an 
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atmosphere, 400 km amagats deep, overlying a thick cloud would explain 
the Stratoscope n pictures. Although the Stratoscope II imagery shows no 
structural detail on the disk, the lack may simply be the result of low contrast 
introduced by the exceptionally wide band-width. 

At infrared wavelengths, Westphal (1972) obtained two photoelectric scans 
across the Uranus disk in the passbands 8 OOOI - 8240^ and 8720l - 8960l . 

Both were obtained simultaneously on 1971 March 8. The atmospheric seeing 
was ~0.5 arc, and a circular 1 arc aperture was used. One scan (8720A - 
8960l) has since been reproduced in the literature {Belton and Vesculus, 1975), 
'but the other remains unpublished. The long wavelength scan is limb-brightened; 
the short wavelength scan shows a flat central region ~ 2 " arc wide together 
with Her D-darkening. Also, the 8340l scan is asymmetric; the west limb 
appears brighter than the east limb. Investigation of the importance of smearing 
effects, due to the natural aperture response of the instrument and atmospheric 
seeing, showed that the tiny (4”arc) disk of Uranus may appear limb-darkened 
even if weak limb-brightening is present. Taking spatial smearing into 
account, Belton and Price (1973) found no obvious contradictions between their 
predictions and the Westphal observations. 

Infrared limb -brightening was confirmed by Sinton (1972). Several 
photographs of Uranus were obtained with a Varo-tube and an 8870l inter- 
ference filter (140l half -width; 220l tenth width) under good seeing 
conditions (~ l".l arc) . Both limb-brightening and polar brightening were 
found, Sinton interpreted the polar brightening in terms of a haze in the 
upper atmosphere, in addition to Rayleigh scattering. Unfortunately, 
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knowledge of the point spread function was too imprecise to permit 
detailed quantitative comparison with the Belton and Price (1973) 
predictions . 

In their critical review of the visual, photoelectric, and photo- 
graphic data, Belton and Vesculus (1975) concluded that markings do 
frequently occur on Uranus. But knowledge of the optical appearance of Uranus 
Uranus is still extraordinarily vague. Thoroughly defining the limitations of 
the semi-infinite, clear H 2 - CH^ atmospheric model would be a significant 
step forward. Results of a systematic attempt to investigate in detail the 
validity of the Belton and Price (1973) predictions are reported in this paper. 

2. OBSERVATIONS 

Uranus is difficult to observe. Being located so far from the sun, the 
planet exhibits only a faint (m^ -i- 6), tiny (~ 4 arc diameter), disk. 

It 

Atmospheric seeing broadening, tjjpically 1 - 2 arc, severely distorts 
the telescopic image. Detecting limb- brightening on Uranus is an extremely 
challenging problem. Fortuitously, photometry of close double stars encounters 
similar observational difficulties which, to a large extent, can be overcome by use 
of a photoelectric area-scanning technique. Pioneering development was 
carried out fay Rakos (1965). initial results were so promising that a scanner 
designed specifically for the simultaneous astrometric and photometric 
study of visual binaries was constructed (Franz, 1966; 1970). With minor 
modifications, the equipment was found to be eminently suitable for the 

Uranus task when mounted at the Cassegrain focus of the 72-inch aperture 
Perkins telescope at Lowell Observatory. 
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Exploring the nature of the Uranus disk profile as a function of geo- 
metrical albedo required the selection of eight narrow waveband (~100A) 
filters spread throughout the region \5500l - \7600l (Fig. 1). Charac- 
teristics of the filters are listed in Table I. Experience showed that choosing 
significantly narrower band widths would lead to inadequate signal/noise ratios 
at the longer wavelengths. To maximize the signai/noise ratio throughout, 
an EMI 9558 (S-20) photomultiplier was employed in the system. 

Selecting the optimum technique for scanning across the Uranus disk 
required careful consideration. Use of a pinhole aperture is one option. 

With excellent seeing, the sensitivity to ifmb- brightening is maximized. 

But two significant problems appear under typical observing conditions. 

First, it is difficult to insure that a pinhole moves precisely across a 
diameter of the Uranus disk. If not, asymmetries in the scan may result. 
Second, theoretical interpretation of the data is extremely difficult even if 
the point spread function is reliably known. Determining the instrumental 
response for a finite circular aperture passing across a circular disk which 
contains a radial intensity gradient is a non -trivial exercise. By comparison, 
slit scans greatly simplify both the observational procedure and the data 
interpretation. If the slit length is significantly greater than the image 
diameter, radially symmetric scans are always assured. It is no longer 
essential for the aperture center to pass precisely through the image center. 
Specifically, accurate determination of the point spread function is readily 
achieved. Moreover, chord integration with slit broadening for either Uranus 
or a stellar image is a tractible analytical problem. 
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TABLE I 
.FILTER SET 


Filter No. 

1 - 

1 Central Wavelength (1) ' 

i' 

1 

5600 

i 2 

i 

6000 

3 

6200 

4 

6400 

5 

6800 

■ 

6 1 

7000 

7 

7300 

1 « i 

7500 


Notes: 1, Filters No. 1 - No. 6 inclusive were manufactured 
by the Optical Coating Laboratory. Tolerances in 
their central wavelengths are +30l. In each case, 
the half -power band width (HPBW) is 100^^201. 

Peak transmission is a nominal 50%. Typical trans- 
mission figures are, 10% at 1.25 HPBW; 1% at 1.8 
HPBW; 0.1% at 3,0 HPBW; <0.1% at > 3.0 HPBW. 

— ■ 2. Filters No. 7 and No. 8 were manufactured by In- 

frared Industries, Thin Film Products Division. 
Tolerances in their central wavelengths are + 15l. 
In both cases, the half -power band width (HPBW) is 
100+151. Peak transmission is a nominal 50%. 


O.ITY OF THE 
■ pnnn 
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Both pinhole and slit scans were used in the observational 

program. Characteristic widths of the pinhole and slit were both chosen 

equal to lOOp (oV645 arc). Wider apertures caused inadequate spatial 
* ■ 

resolution. Narrower apertures gave a poor signal/noise ratio throughout 
the wavelength range under study. Scan orientations were either north- 
south or east-west. NS scans were preferred; image broadening, through 
telescope drive error, could not occur. For EW scans, individual filter 
integration times were kept correspondingly short. Table n gives the 
complete Uranus observing log. 

3. RESULTS 

3. L Coarse Analysis 

Expeditious processing of the Uranus data required a simple, objective, 
method for detecting the coarsest features of the true intensity profile of the 
disk. Whenever seeing broadening is significant, limb-brightening on Uranus 
will manifest itself by a general widening of the planetary image. Similarly, 
limb- darkening will cause a relative image-narrowing. All observational 
data were acquired under ordinary Flagstaff atmospheric seeing conditions, 

rt 

stellar image diameters (half -power) being typically ~ 2 arc. Both pinhole- 
and sUt-scans of Uranus in all wavebands showed a characteristic Gaussian 
shape. Determining the Gaussian 1/e-width of the Uranus image in each 
waveband should therefore provide a reliable method for comparing relative 
limb-darkening/Iimb-brightening effects . 

Coarse analysis of the observations consisted of determining, by the 
least squares technique, the Gaussian 1/e-width for each and every integrated 


TABLE n 
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URANUS OBSERVRTG LOG 


1 

SCAN 

1 

SKY 

i 

: 

DATE(U.T.) 

TYPE 

FILTERS ! 

TRANSPARENCY 

SEEING 

REMARKS 

1974 May 1 

‘P 

■ ■ 1 

1,2, 3,4, 5, 6, 8 1 

5 

2-3 

E<-»w ; 

1974 May 26 

P ■ 

1,2, 3,4 ! 

5 

3-4 

N«-^S ; 

1974 June 5 

P 

2,3,4 

4 

2 

N^S ; 

i 

1975 May 24 ' 

S 

1,2,3, 4, 5, 6,7, 8 

4-5 

’ 

3-4 

N<-^S; Scale * 
known but no 
PSF informa- 
tion 1 

j 

1975 May 27 

S 

1,2, 3,4,5, 6, 7, 8 

5-3 

' 

3-2'" 

I 

N<-^S;E^W; 
Scale from 
ADS 9053;PSF 
from<^ Virginis 

; 1975 June 15 

s 

' 

' 

1,2, 3, 4, 5, 6, 7, 8 

j 

5 

1 

1 

3 

E«^W;Scale 
known; PSF ; 

from K & i j 

Virginis 1 

( 

1975 June 16 

s 

1 

i 

1,2,3, 4, 5, 6, 7, 8 

5 

3 

E^WjScale ' 
known; PSF 
fromtVir- J 

ginis 1 

1975 June 16 

i 

P 

6,7,8 

5 

3 

E^W 

1975 June 17 

P 

1,2, 3,4, 6, 7, 8 

5 

3 

N^S 

1 

j 1975 June 19 

I 


2,3,7, 8 

5 

i 3 

1 

E^W 


Notes: l. Scan type is either pinhole (P) or slit (S) 

2, Filter No.'s are taken from Table t 

Sky transparency and seeing conditions are each given on a scale 0-5 
(i.e. Worst-Best) 


4. Information concerning scale calibration and seeing data are given in 
the Remarks column. 
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pinhole- and slit-scan. Subsequent numerical experiments showed that 
the use of more sophisticated profile shapes had a negligible effect on 
the results Franz et al (1971) and Franz (1973). Data secured using 
different scan-types and scan-orientations were initially kept segregated 
to determine if the results were sensitive to the mode of observation. 
Specifically, data taken on a particular night, by a particular scan technique, 
were treated as a complete set. Each set was obtained by a rapid, symmetrical, 
forward-backward cycling of the filter wheel. Integration times for each filter 
were normally in the range 20 - 500 seconds. For all observations within 
a given set, the seeing and transparency could be taken as constant. Detailed 
point spread function information, and accurate linear scale calibrations, 
were not always available on each night. Consequently, only relative profile 
widths were determined within each set; all profiles were conveniently 
normalized to the mean for Filter No. 2 (XOOOOA). Relative Gaussian profile 
widths are listed in Table HI. The rms error of the individual relative width 
determinations are also shown. Where no rms error is available due to 
insufficient data, the uncertainty is indicated by a question mark. In 
parentheses, following each rms error, is the number of individual width 
determinations from which the mean value was calculated. 

Table IV summarizes the results obtained from each combination of 
scan -type and scan-orientation. Corresponding mean values for the relative 
Gaussian widths have been derived from the values Listed in Table HI. 

Individual rms errors quoted in Table IV are based on the spread of the 
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TABLE m 

URANUS DISK: RELATIVE GAUSSIAN PROMLE WIDTHS 


"'-"'-^-.^LTER NO. 
DATE (U. 

1 

2 

3 

4 ! 

i 

5 

6 

7 

8 

Pinhole Scans(N^S) 




1 





1974 May 26 i 

1.01+? (1) 

1.00+.02(3) 

1. 06+. 02(2) 

1.00+? (1)| 

X 

X 

X 

X 

1974 June 5* 

X 

1.00+.00(3) 

1. 02+. 01(3) 

0. 99+. 01(3) 

X 

X 

X 

X 


X 

1.00+. 02(135 

1.04+.030.4) 

0. 98+. 0302(1 

X 

X 

X 

X 

1975 June 17 

0.99+? (1) 

1.00+.03(4) 

1.03 + .03(5) 

1.04 +.03 (2) 

X 

1.01+? (1) 

1.15+.02(2) 

1.01+? (1] 

Pinliole Scans(E^W) 


1 

1 







1974 May 1'*' 

0.98 +.02 (3) 

1.00j-.03(6) 

1,03^.03(7) 

1.02 + .04(5) 

1.04+.03(3) 

1.05+.01(3) 

X 

1.04+.0K2) 


1.01+? (1) 

1.00 +.01 (4) 

i.04+.02(8) 

1.00 +.02 (4) 

1.00+? (l)i 

0.99+? (1) 

X 

X 

1975 June 16 

X 

X 

X 

X 

X i 

1. GO +.02 (2) 

1.28+? (1) 

1.01+ ? (1) 

1975 June 19 

X 

l,00+.06(6) 

0.99+.06(4) 

X 

X 

X 

1.06+? (1) 

0.89+? (1) 

Slit Scans (N<-^S) 







1 


1975 May 24 

Il.02+.Q1(2) 

1.00 +.00(2) 

i 1.03 + .00(2) 

0.98+.01(2) 

1.00 +.00 (2) 

0.99+.01(3) 

1.08 +.02 (4) 

0,98+.01(2) 

19(0 May 27 

0.99 +.01 (2) 

1.00+.00(2) 

1.02 + .01{2) 

0.99+.01{2) 

0.96+.01(2) 

0.95+.00{2) 

1.05+.01(2) 

0.97+,01(2) 

Slit Scans (E<i-^W) 









1975 May 27 

0.96 + .00 (2) 

1.00+.01(2) 

1.00 + .03(2) 

0.97+.02(2) 

0.95+.03(2) 

0.97+.03(2) 

1.07+.01(2) 

0.95+.02(2) 

1975 June 15 

1.04+.03(3) 

1.00 +.06 (4) 

U03+.03{4) 

1.02+.04(5) 

1.02+.06(4) 

1.01+.03{4) 

1.08 +.03 (4) 

i.01+.01(5) 

1975 June 16 

0.99+.OK2) 

1.00+.01{2) 

1.02±.01{2) 

1,04 +.03 (2) 

, 1.02+.02(2) 

1.03+.02{4) 

1.08+.01(2) 

0.99+.01(2) 


Notea: * Scale change between two sets of measui’ements 

f Cam change between two sets of measurements t? 


TABLE IV 


URANUS DISK: GAUSSIAN WIDTHS' SUMMARY 


^^JF-ILTER NO. 

date(u.t7t\^ 

1 

2 

3 

4 

5 

6 

7 . 

8 

Pinhole (N<-^S) 

l.OO+.Ol 

1,00 

1.04+.02 

1 . 00 f .02 

X 

1 * 01 + ? 

1.15+.01 

1 . 01 +? 

Pinhole (E^W) 

1 . 00 J -.02 

1.00 

1 . 02 +.02 

l.Ol+.Ol 

1 . 02+.02 

1.01+.03 

1.17+.11 

0.98+. 06 

Slit (N<-^S) 

1 , 01+. 02 

1.00 

1.03+.01 

0.99+.01 

0.98+.02 

0.97+.02 

1.07+.02 

0.98+.01 

Slit (E<~^W) 

1.00+.03 

1.00 

1 . 02+.01 

1.01+.03 

1.00+.03 

1.00+.03 

1.08+.01 

0.98+.03 

Mean 

l.OO+.Ol 

1.00 

1.03+.01 

l.OO+.Ol 

l.OO+.Ol 

1 . 00+.02 

1.12 +.04 

0.99+.01 


Notes; 1. Each type of scan data is normalized to Filter No. 2 

2. If no rma error is available through insufficient data, the 
uncertainty is indicated by a question mark. 
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relevant values listed in Table m. Throughout, ail mean relative 
Gaussian widths have been normalized to Filter No. 2. For each filter, 
a final relative Gaussian width has also been derived by averaging the 
values listed in Table IV itself. The corresponding rms error has been 
obtained from the actual spread in the Table IV numbers. Several con- 
clusions can be drawn: 

1. Except for Filters No. 3 and No. 7, there is no evidence of 
any change in the Gaussian 1/ e width of the Uranus image with 
wavelength. 

2. Filters No. 3 and No. 7 both show evidence of an increase in the 
Gaussian 1/e width relative to Filter No. 2. Moreover, the 
increase is larger for Filter No. 7 than for Filter No. 3. 

3. For Filter No, 7, the relative Gaussian 1/e width is significantly 
larger if pinhole- scanning rather than slit-scanning is employed. 

4. Scan-orientation does not signiihcantly affect the relative 
Gaussian l/e widths. Evidently, the telescope guiding was 
accurate. 

Our results are most readily explained by the presence of Limb-brightening 
in the wavebands corresponding to Filters No. 3 and No. 7. Moreover, the 
limb-brightening should be more pronounced for Filter No. 7 than for Filter 
No. 3. 

Composite pinhole NS scans were constructed to maximize the visibility 
of limb -brightening on the Uranus disk. Specimen results are illustrated 
in Fig. 2 . Observations taken on 1975 June 17 were selected on the basis 
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of stability both in the atmospheric seeing quality and in the sky trans- 
parency^ Filter No, 7 corresponds to the deep X7300I. CH^ band in which 
limb -brightening is to be expected. Filters No. 6 and No. 8 refer to 
adjacent wavebands in which limb-brightening is not anticipated. The Filter 
No. 7 profile was constructed by summing two integrated scans made up of " 
1000 and 500 individual one-second scans respectively. Integrated scans, 
each made up of 200 individual scans were used to form the composite of filters 
No. 6 and No, 8. Gaussian curve-fitting parameters were used in the co- 
location, background elimination, and normalizatic.n, of the integrated scans. 
For Uranus, the photodetector response for filters No. 6 and No, 8 was 
6 times greater than for filter No. 7. Significant differences in the signal/ 
noise ratio are apparent in Fig. 2. A difference in the profile widths can 
be detected through inspection; Gaussian curve-fitting shows the filter No. 7 
profile to be significantly wider Table in. 

3.2. Fine Analysis 

Detailed quantitative studies of selected Uranus observational data 
were carried out to confirm and enhance Che results of the coarse analysis. 
Attention was concentrated on the X7300l CH^ band where, limb-brightening 
appeared to be most pronounced. Determining whether the apparent 
limb -brightening was in fact absolute, or merely relative, was of prime 
interest. Results obtained would have a direct bearing on the interpretation 
of the X6190l CH^ band data. Several criteria were used in the data selection. 

First, stability in the observing conditions was of paramount importance. 
Not only must the sky transparency be uniformly excellent, to maximize 
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photometric accuracy, but the atmospheric seeing must be exceptionally 
stable. While excellent seeing would have been ideal, constancy raiher 
than narrowness was the watchword for the point spread function. Second, 
accurate scale-calibration of the linear scan was essential. Third, slit- 
scans only could be utilized; quantitative interpretation of pinhole data was 
intractible. Fourth, since maximizing the Umb-brightening/limb-darkening 
contrast was vital, composite profiles for Filter No. 7 and for Filters 
No. 6 and No. 8 were considered to be essential for detailed comparison. 

Their construction required the availability of significant numbers of inte- 
grated scans. Observations made on 1975 June 15 were found to be optimum 
_cf : Tables n and El. Scan -orientation did not compromise the selection. 

NS- or EW-scans had been found to lead to essentially identical profiles. 

Point spread function. (PSF) data, in the relevant wavebands, were obtained 
from slit-scans of two stars, f< Vir and L Vir, both near Uranus. 

To validate atmospheric seeing stability, the Uranus and stellar scans were 
intermingled. Eleven integrated stellar scans, taken through Filters No. 6 
(3), No. 7 (4), and No. 8 (4), were used to derive a composite PSF profile 
relevant to the quantitative interpretation of the Uranus data (Fig. 3). No 
significant color-dependent differences in the shapes and widths of the stellar 
profiles were apparent. Each integrated stellar scan was uniformly comprised 
of 20 one-second sweeps. Gaussian curve-fitting parameters were used in the 
colocation, background elimination, mirror-imaging, and normalization, 
of the stellar data. The linear scale was determined using the visual 


■binary 5 Ser (ADS 9701: Separation {1975,5) = 3,96_+0.02 arc). Quantita- 
tively, the atmospheric seeing stability can be described by the rms 
deviation of the eleven individual Gaussian l/e-width parameters of the 
chosen stellar scans. An rms error of only 5.8 percent was derived. 

Selected Uranus data comprised 13 integrated scans spread throughout 
the three wavebands in the following manner; Filter No. 6 (4), Filter 
No. 7 (4), and Filter No, 8(5)., Each Filter No.. 6 scan consisted of 20 
one-second sweeps. Three Filter No. 7 scans each consisted of 200 sweeps; 
, a fourth scan was made up of 100 sweeps. Each Filter No. 8 scan consisted 
of 100 sweeps. Derived Gaussian l/e-widths for all Filter No. 6 and Filter 
No. 8 integrated scans were statistically intercompared to determine if the 
selected waveband, or the number of individual sweeps within each scan, 
significantly affected profile shape. No indications of any such differences 
were found. A composite "continuum" profile was therefore derived by 
taking all Filter No. 6 and Filter No. 8 integrated scans together. For 
FiEer No. 7, the number of sweeps making up each integrated scan did not 
affect the Gaussian l/e-width. A composite \7300A CH^ band profile was 
therefore derived by taking all Filter No . 7 integrated scans together . 
Expeditious use was made of the Gaussian curve-fitting parameters to 
derive each composite Uranus profile. Special care was taken to insure 
that each composite scan was normalized at the precise center of the disk. 
Final normalization utilized data within 5 channels (0 Jl5 arc) radius of the 


disk center . 
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The composite Uranus sUt^scan profiles in the X7300i CH^ band, and 
in adjacent "continuum" regions, are compared in Figs. 4 and 5. Brief 
inspection of the observational data shows the CH^ band profile to be 
significantly broader. Relative broadening is particularly noticable at 
distances 1-2 Uranus angular radii from the disk center. This is precisely 
where effects of limb-brightening would be expected to show to maximum 
advantage. Relative limb-brightening, first indicated by the coarse analysis, 
is confirmed. 

Inferring the true intensity profiles of the Uranus disk can be based on 
either of two techniques. One may attempt to deconvolve the various smear- 
ing functions affecting each composite scan. Alternatively, numerical 
prediction of the observations, based on subjecting individual models of the 
true disk profile to various smearing mechanisms, may be used. For 
convenience, the numerical prediction technique was adopted. Three distinct 
models of the Uranus disk were considered- Model I represents a disk of 
uniform surface brightness. Model II assumes limb-brightening; the surface 
brightness varies linearly from zero at the disk center to a finite value at 
the limb. Model ni assumes limb-darkening; the surface brightness varies 
linearly from a finite value at the disk center to zero at the limb. In modelling 
the Uranus disk, circular symmetry was adopted. From the results of 
Danielson, Tomasko, and Savage (1972), the Uranus radius was taken as 
25,900 kms. Distance from the Earth to Uranus at the time of observation 
was taken from the 1975 American Ephemeris and Nautical Almanac to be 
17.89 A.U. The corresponding angular diameter of Uranus (unbroadened) 
was therefore 3.99"' arc. 
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Three distinct smearing mechanisms together affect the original 
Uranus disk to produce the observed planetary profile. First, every point 
on the disk is subjected to two-dimensional seeing broadening by a circularly 
symmetric .point spread function. Second, chord- integration of the radially 
symmetric planetary image occurs because of the type of scan employed. 
Third, slit-broadening of the chord -integrated profile results from the need 
to use a finite scanning aperture. Chord-integration and slit-broadening both 
affect the original PSF shape to produce the observed slit scan of a stellar 
image. For convenience, the original PSF function was taken as Gaussian. 

’ Not only does a radially symmetric Gaussian profile provide an adequate 
first order description of the original point spread function, but its adoption 
leads to important mathematical simplifications in the theoretical convolution 
of the Uranus disk. 

PSF observational data and theoretical predictions are compared 
in Fig. 3 . Observational points refer to the composite slit-scan profile 
derived from individual stellar images. Theoretical predictions utilize 
the known slit width (0V645 arc) together with individual selected values 
for the Gaussian 1/e-width (o) of the original PSF . In keeping 
with the experiences of Franz et al (1971) and Franz (1973), no unique 
Gaussian curve provides a perfect fit to the observations. Nevertheless, 
limiting values for an effective Gaussian 1/e-width can be derived. By 
inspection of Fig. 3 , it must lie in the range 1V25 s a ^ 1V75 . Evidently, 
a narrower Gaussian curve is required in the core of the stellar image than 
in the wings. The optimum Gaussian 1/e-width (ct) was taken as I’lS arc. 
Extreme error bars -0V25 were adopted. 
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Determining the reality of absolute limb-brightening in the CH^ band 
■was the objective for the theoretical predictions shown in Fig. 4, Model I, 
a uniform circular disk, was selected to describe the unbroadened Uranus 
image. Detailed image convolution was carried out using three 
individual choices for the Gaussian l/e- width (a) of the original PSF i.e. 

1 V 25 , iVs, and 1 V 75 . In Fig. 4, the central curi'e refers to the intermediate 
Gaussian l/e-width (iJs). CE^ band observations lie significantly above 
the central curve, strongly indicating that the apparent limb-brightening is 
in fact absolute. Even if the widest permissible Gaussian l/e-width is 
adopted (top curve), the Uranus disk must still be uniformly bright. If the 
narrowest PSF is adopted (lower curve), the need for absolute limb- 
brightening becomes very obvious. For the so-called continuum profile, 
use of the intermediate Gaussian l/e-width strongly suggests that the 
unbroadened Uranus disk is uniform. Use of the widest PSF indicates 
limb -darkening; the narrowest PSF leads to limb -brightening. 

Determining the sensitivity of the slit-scan profile to extreme variations 
in the original distribution of surface brightness on the Uranus disk was the 
objective of the theoretical predictions shown in Fig. 5 . Three distinct 
intensity profiles were selected to describe the Uranus disk -viz. Model I 
(Uniform disk), Model n (Extreme limb-brightening), and Model DI (Extreme 
limb-darkening). In each case, the original PSF was defined by the intermediate 
Gaussian l/e-width (1.5). In Fig. 5, Models I, II and HE are represented by 
the middle, top, and lower curves respectively. Evidently, extreme limb- 
brightening must be present on Uranus in the X7300l CH^ band. If not, the 
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phenomenon would never have become apparent under the mediocre seeing 

d du2!*ixi*^ the course of the observations. Under 
better seeing conditions, limb-brightening should become very striidng. 

4. CONCLUSIONS 

From our coarse analysis, the Uranus disk must exhibit limb-brightening 
in two CH^ bands (\6190A; \ 7300 A) relative to all other wavebands studied. 
Whether the apparent limb-brightening was in fact absolute, or only relative, 
could not be immediately determined. Except for the two CH^ bands, the 
■ shapes and widths of the Uranus images were essentially identical. Since 
limb -darkening is known to occur at visual wavelengths (X5600A), it must 
occur in all so-called continuum wavebands studied. No systematic wave- 
length-dependent effects were found in the original point spread function which 
might have affected interwaveband comparisons, PSF data for all wavebands 
showed no significant variation in the representative Gaussia:. 1/e-width 
with wavelength. 

From our fine analysis, absolute limb-brightening was found in the 
X7300l CH^ band. Taking into account results of the coarse analysis, absolute 
limb-darkening was found in the adjacent wavebands. No definite conclusions 
could be drawn regarding the absolute nature of limb- brightening 
in the X6190^ CH^ band. If absolute limb-brightening does occur, it must 
be considerably less pronounced then in the X7300A CH^ band. 

Although our results are consistent with the basic predictions of Belton 
and Price (1973), the validity of the H 2 - CH^, semi-infinite, clear atmosphere 
remains to be proven. First, further observations secured under improved 
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seeing conditions are required to better define the Uranus disk profiles. 
Second, alternative explanations of the limb-brightening phenomenon 
must be thoroughly evaluated. For instance, the observational consequences 


of a thin aerosol haze, located high in an otherwise clear - CH. 
atmosphere, require detailed prediction. Effects introduced by the 
presence of a deep dense NH^ cloud layer must also be considered. Further 
observations and theoretical analyses are under way. 
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Fig. 1. The schematic variation of the Uranus geometrical albedo, 
p, with wavelength, X (after Younkin, 1970). The values of geometrical 
albedo given by Younkin have been adjusted to the radius derived by 
Danielson et al (1972). Limb-darkening is predicted to occur in the shaded 
regions; limb-brightening in the unshaded regions (Belton and Price, 1973). 
The multi-color filters employed in the observations are shown by the 
numbered bars. Approximate extremities in their transmission curves are 
indicated. 

Fig. 2. Composite NS pinhole scans of Uranus obtained on 1975 June 17. 
Each point corresponds to a discrete channel of the linear scan. 

Fig. 3. The Composite PSF slit-scan profile . Theoretical curves are 
based on individual choices for the original PSF Gaussian 1/e-width (o). 

Slit width is 0.645 arc. 

Fig, 4. The composite Uranus slit-scan profiles. Filter No. 7 (CH^ 
band) and Filter No. 6/No. 8 (Continuum) data are plotted separately. 
Theoretical predictions refer to a uniform planetary disk. Individual choices 
for the original PSF Gaussian 1/e- width (cr) were 1V25, 1V5 and 1V75. The 

slit-width was 0V645 arc. _ 

Fig. 5 . The composite Uranus slit-scan profiles. Filter No. 7 

(CH. band) and Filter No. 6/No. 8 (Continuum) data are plotted separately. 
Theoretical predictions refer to the intermediate original PSF Gaussian 
1/e-width (175) coupled with individual choices for the planetary disk 
profile vix. Model I (Uniform Disk), Model n (Extreme Umb-brightening) 
and Model in (Extreme limb-darkening) . Slit-width was 0V645 arc. 







REUVTIVE INTENSITY 









RADII 



{"ARC) 





RELATIVE INTENSITY 







APPENDIX E 


MARS: 

TOPOGRAPHIC CONTROL OF CLOUDS, 
1907-1973 


by 


William K. Hartmann 
Planetary Science Institute 
2030 E. Speedway, Suite 201 
Tucson, Arizona 85719 


Receiyed 


Revised 


No. of Copies: 4 

No. of MS Pages: 19 

No. of Figures: 2 


No. of Tables: 


0 



Proposed Running Head; 


Proofs should be sent to; 

Dr. William K. Hartmann 
Planetary Science Institute 
2030 East Speedway, Suite 201 
Tucson, Ariaora 85719 


E-3 


ABSTRACT 

Mariner 9 high-resolution photos and topographic information 
have been used to analyze '’blue" and "red” clouds reported over many 
decades. Photographic records of cloud positions, in the Lowell 
Collection, covering a 66-year period from 1907 to 1973, are compared 
with topographic features. 

A sample of 77 "blue" cloud sites lay preferentially at the highest 
Martian elevations; 60% centered precisely on the 7 major volcanic 
mountain peaks; another 16% lay on substantial slopes or contacts be- 
tween cratered terrain and lower plains. Median altitude of blue cloud 
sites was 2. 1 km above the global topographic median. These results 
support the conclusion that most blue clouds are orographic uplift clouds, 
composed of condensates. 

Over half of 131 sporadic yellowish or red clouds, identified with 
filter photography, were associated with blue clouds or volcanoes, and 
thus probably did not represent dust storm phenomena, contrary to a 
commonly held belief. Of 88 "possible dust clouds," about 2/3 occur 
at borders between light and dark areas, in the light regions, perhaps 
sites of current depositional or denudational activity. Median altitude 
of "possible dust cloud” sites was 0 . 5 km below the global topographic 
median. Major dust storms begin in a few core areas, two of which 
associate with major basins of Hellas and Isidis. Basins are probably 
reservoirs of mobile dust. 


/s 
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I. INTRODUCTION 

Through the years, many researchers have tabulated the behavior 
of clouds and other variable features on Mars (Martin and Baum, 1969; 
Capen, 1974; Martin, 1974) and even attempted to infer physical 
processes and topographical characteristics from such studies (Wells, 
1966; Arvidson, 1972; Schlosser and Haupt, 1973). However, . only 
with the combination of synthesised global topographic information — 
not available in relatively final form until 1975 — and high-resolution 
Mariner 9 global photography, could variable feature behavior be 
reliably correlated with Martian topography and geology. 

The purpose of this study was to seek systematics of Martian 
clouds as a function of topography and probable cause. The study was 
particularly designed to utilize archives of observations that stretch 
back nearly 70 years, thus providing a more general survey than the 
near-hourly photography of recent phenomena, made possible by the 
world-wide planetary patrol photography of recent years (Capen, 1974; 
Martin, 1974). This is believed to be the first such long-term survey 
with reliable global topographic information. 

Much effort was spent over the last 70 years to understand these 
observed phenomena. Our interest in modern, vastly more detailed 
spacecraft observations should not make us forget our opportunity to 
"clean up loose ends" by determining what phenomena the older earth- 


baaed observers were actually seeing. 
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E. MARTIAN TOPOGRAPHY 

During this study, a major synthesis of topographic information 
(occultation, radar, pressure surfaces) was being completed by the 
U.S. Geological Survey. Two versions, August 1974 and June 1975, 
were consulted, and the conclusions of this study were not affected by 
the modest changes between these versions in the non-polar regions 
, studied here . The final tabulations and analyses of data were made 
with the T.25,000,000 topographic map (M 25M 3RMC; June 1975). 

It includes topographic contours, structural features in shaded relief, 
and albedo markings taken from earth-based and Mariner 9 imagery. 
The reference elevation (zero km) is defined as the 6,1 mb pressure 
level. Although topographic data continue to be refined, the major 
conclusions of this paper are likely to remain the same. 

Elevation distributions on Mars were studied in order to compare 
with the distributions of variable features such as blue clouds and red 
clouds. M^rs was divided into 720 equal-area blocks, from which 
mean elevations and a histogram of elevation frequency were obtained, 
following the technique used in an earlier study (Hartmann, 1973). 
Figure 1 compares this histogram, or hypsometric diagram, with the 
elevation distribution of blue clouds and red clouds, discussed below. 

A median elevation for the entire surface of Mars was found to 
be +2.4 km; this will be compared with median elevations for clouds. 
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m. BLUE CLOUDS 

The primary source of information for this study was an unpub- 
lished Lowell Observatory report by Martin and Baum (1969). The 
authors reported a search of 5,000 Mars plates or any groups of plates 
on which discrete clouds were observed over a period of a few days. 

The purpose of that search was to measure rates of motion in clouds, 
and so there was a selectional bias toward discrete clouds whose motions 
could be traced. This tended to restrict the study to well-defined bright 
clouds, a selection which appears useful for the present study. Martin 
and Baum reported 95 cloud histories ranging from 1907 through 1958. 
The vast majority of these clouds were photographed on blue -sensitive 
or yellow-green-sensitive plates and are believed to be blue clouds or 
bluish-white clouds. 

Martin and Baum described cloud histories against a background 
of traditional pre-Mariner maps of albedo features, and some interpreta- 
tion was required by them, and by me in a few cases. For example, 
in May 1937 a bright cloud was seen in successively northeastern 
positions in apparent motion across the Maria Tyrrhenum region. Be- 
tween the fourth and fifth day the cloud appeared to shift some 240 km 
to a new position over Elysium. Martin and Baum noted that this might 
imply a new Elysium cloud rather than a continuation of the earlier 
cloud, but listed the whole group as a single cloud history. Noting that 
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the final cloud position precisely overlapped the Elysium volcano, I 
chose to plot it as a separate cloud position. 

I also deleted 11 clouds observed primarily on red-sensitive 
plates and several other clouds of ambiguous color, observed with 
different filters. In this way the Martin-Baum list was altered to 77 
blue cloud positions. Positions plotted were essentially the initial 
, positions or area of origin of the blue cloud, as the purpose was to 
study geographic conditions at sites where clouds originate. 

Of the 77 clouds, 46 or 6Q% formed directly over the flanks or 
summits of 7 major volcanic mountains on Mars . These were: Olympus 
Mons, Ascraeus Mons, Pavonis Mons, Arsia Mods, Tharsis Tholus, 

Alba Mons, and Elysium Mons. In addition, another 1Q% occurred on 
noticeable slopes (defined by noticeably higher than average contour 
line density — tjrpically slopes averaging ^ 2 km rise over 500 km, or 
2: 0.4%), contacts between higher cratered areas, and lower featureless 
plains, or basin rims. 

The very strong correlation with 7 volcanic mountains standing 
about 2 to 16 km above their surroundings, whose existences were un- 
known at the time the blue clouds were mapped by Martin and Baum, 
strongly supports the explanation of these clouds as products of con- 
densation of water ice due to orographic uplift or possible local degassing. 
The latter interpretation was developed for similar clouds in the same 


positions, observed independently by Mariner 9 (Leovy, Briggs, and 
Smith, 1973), and the cloud composition was determined spectro- 
scopically to be water ice. The importance of orographic uplift for 
blue clouds is similarly supported by the association of the additional 
16% with slopes and contacts between structural units at different 
altitudes. Thus 3/4 of all blue clouds can be described as associated 
. with volcanoes or sloping contacts, while the other 1/4 showed no clear 
systematic pattern. 

The median elevation of the sites of 77 blue clouds was found to 
be 4.5 km, significantly higher than the global median of 2,4 km. 

This again indicates the association of blue clouds with high-altitude 
air masses, and confirms the conclusion of Schlosser and Haupt (1973) 
who found that evening and morning terminator clouds correlate with 
high radar-determined elevations in the equatorial regions. 

IV. possible dust clouds 

To clarify the rest of the study, cloud terminology should be 
more clearly defined , Many photos in the Lowell Collection are made 
with specific filters, usually blue, yellow, or, less commonly, red. 
Clouds seen in such photos have been called "blue,” "yellow," or "red," 
giving rise to a mistaken impression that blue and red clouds are dis- 
tinctly different. For example, while blue clouds are often associated 
with condensates, red clouds are often regarded as dust clouds. 
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In reality, the types are mixed. The tabulation fay Martin and 
Baum (1969) and the original plots of clouds from 1907 to 1958, pre- 
pared fay Martin for that report, show that clouds detected in yellow 
light are frequently found to fae identical with clouds shown in blue 
light, if blue plates are available for the date in question. Using 
Martin's original plots (kindly provided by Martin for this study), I 
identified and mapped 131 clouds detected in yellow or red light, 
plotting initial positions when motion was detected. Of these, 52 (40%) 
could be identified with major blue clouds or blue cloud complexes. 

An additional 18 (14%) lay on summits of major volcanoes, and were 
assumed to be identical with blue clouds, in the absence of other obser- 
vations . To obtain a sample that might more likely be identified with 
dust clouds, these 70 clouds (54% of the sample) were rejected, leaving 
a sample of 61 yellow or red clouds that will be called "possible dust 
clouds." 

No prominent simple relation was seen between the "possible dust 
cloud" sites and Martian topography or structure, although 21 clouds 
(34%) overlapped the floors and rims of basins or the Valles Marineris 
complex , 

The elevation distribution of the 61 "possible dust clouds" is shown 
in Figure 1. The median elevation is +1.9 km, falling 1/2 km below the 
globe-wide median of 2,4 km. The lower than average elevation of 



E-10 


^'possible dust clouds'* is consistent with experimental findings of 
Sagan and Bagnold (1975) who concluded that raising of dust particles 
(especially 1-100 micrometer diameter) by winds is easiest in low, 
cold regions . 

Topographic correlation is not as striking for these red and 
yellow "possible dust clouds" as for the blue clouds, however. 

. Figure 2 shows the geographical distribution of the 61 "possibV dust 
clouds," augmented with 27 other red clouds taken from Sl.ipher (1962, 

12 clouds from 1922 to 1956 not included by Martin and Baum, 1969), 

Capen (1974, 1 cloud), and IVIartin (1974, 1976; clouds observed in 
the 1971 and 1973 dust storms). Here clusters are found near the 
margin of Hellas and Isidis basins, but the correlation is not striking. 

Of the 88 clouds about 20% are associated with the Hellas and Isidis 
basins, 

A more interesting correlation is that the "possible dust clouds" 
tend to be located along dark/light albedo boundaries. In the light regions . 
In Figure 2, of the 88 clouds about 66% may be said to meet this criterion; 
about 75% appear to be associated with borders, either on dark or light 
sides, though it is difficult to define these criteria exactly. If these 
88 "possible dust clouds" are really sites of wind-blown dust, this 
finding may indicate that dust is more mobile along the light sides of 
albedo boundary, either because of dust's mechanical properties there 
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or because of wind properties there. (It should be noted again that 
these 88 clouds are defined by being red or yellow, not coinciding with 
blue clouds, and not located on volcanoes.) This finding is consistent 
with the Mariner 9 result that classical albedo markings are associated 
with dust deposition or denudation, as markings change shape from 
year to year under current Martian conditions . 

V. MAJOR DUST STORMS 

To the observations of sporadic dust clouds can be added data on 
dust raised during beginnings of major dust storms. 

Leonard Martin, of Lowell Observatory, has mapped on a daily 
basis the positions of clouds observed during the outbreak of the 1956, 
1971, and 1973 major dust storms. Summaries of much of this work 
have been published (Martin, 1974, 1976). 

L. Martin kindly made available to the author his original plots 
of yellow dust cloud motions on a daily basis for the 1956 storm, and 
on an hourly basis for the 1971 and 1973 storms. Analysis of daily 
positions of clouds on Mars topographic maps, one map for each year, 
showed that certain nuclei exist for the origin and early development 
of Martian dust storms. As the storms spread across the planet, the 
nuclei may become active in advance of the front. 
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The three main nuclei are marked by white triangles in Figure 2, 
which should be distinguished by their greater activity from the black 
spots, or sites of clouds observed sporadically over the last seven 
decades. The two major storm-spawning nuclei in the Hellas and 
Isidis basins are consistent with the idea that storms are favored to 
begin in low basins that act as reservoirs of mobile dust. This idea 
does not account for the third major core area, lying on the high plain 
of Solis Planum, at elevations around +8 km. The nuclei in which dust 
storms are spawned presumably are related both to topography and to 
patterns of high-speed winds. The theoretical dust storm model of 
Gierasch and Goody (1973), for example, predicts that major self- 

I 

sustaining dust storms (as opposed to simple dislodging of dust by 
wind, studied by Sagan and Bagnold, mentioned above) will form "ui 
a high plateau area at times of maximum insolation.*’ The Hellespontus 
rim of Hellas and the Solis Planum region both match the ’’high plateau” 
description (Gierasch and Goody noted the first but not the second.) 
Other workers have calculated wind flow patterns and concluded that 
especially high winds occur on the inner slopes of Hellas basin and 
Isidis basin, consistent with these two regions being storm- spawning 
cores {Sagan, private communication; Mass and Sagan, in press), 

Sagan and Bagnold conclude that Hellas is the most likely place to raise 
dust on Mars, because of the high pressure on its floor, because it has 
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the highest vertical velocity components predicted on the planet. To 
this we may add that Hellas is also a probable major reservoir of 
dust, based on photo-interpretation (e.g. Hartmann, 1973). 

These results can be compared with results of Arvidson (1972) 
who analyzed 58 initial positions of yellow clouds originally tabulated 


from various observations gathered 

statistical clusters in: 


Syrtis Major - 
Isidis Region 

280, +15 

Hellespontus - 
Hellas 

310, -35 

Noachis Region 

340, -40 

Eos Region 

o 

1 — t 
1 

o 

with secondary highs in: 


Mare Sirenum 

140, -40 

Solis 



by Wells (1966). Arvidson finds 


Matches Isidis core, 
Figure 2, '^0 km 

Matches Hellas core, 
Fipire 2, -1 km 

Cratered upland, +4 km 

E end Valles Marineris, 
chaotic, ~0 km 


Cratered plain; 4-4 km 

Matches Solis Planum 
core. Figure 2, -^6 km 


Comparing Arvidson' s study with the present work, one concludes 
that the Hellas and Isidis basins, with their high winds, are major sites 
for producing dust storms, and probably are reservoirs of mobile dust. 



VI. SUMMARY 


Tlie major inferences from this iwork are the following. Classical 
blue clouds are mostly orographic condensation clouds associated with 
mountains and slopes, with 60% lying on Martian volcanoes. Many 
classical yellow clouds are identical with the blue clouds. Of other 
yellow-red clouds, many are true dust clouds generated in low areas, 
t They appear to occur preferentially near boundaries between dark and 
light albedo regions, in light areas. Two interpretations of the latter 
finding are possible. First, it may be an observational selection 
effect, whereby clouds are detected more often if they are raised near 
borders, where they may contrast with coloration of adjacent areas. 
Alternatively, the finding may give information about the physical state 
of the albedo borders, Indicating either excessive winds or high mobility 
of dust in these regions. Origin and spread of major dust storms in- 
volves activity in a few "core" areas, such as Hellas and Isidis basins, 
which may be reservoirs (potential wells) where mobile dust collects. 
Other "core areas" have less obvious topographic distinction and the 
reason for their association with storms is leas apparent . 
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ELEVATION (KM) 




Median 1.9 


Figure 1. Altitude distributions of general Martian surface (top), 
sites of 77 blue clouds (middle), and sites of 61 (red) 
"possible dust clouds" (bottom). Shaded blocks repre- 
sent blue clouds over volcanic mountains. 




All I 


o 


© 




& 


© 


© 


© 


o 


© 


o 


Q 


© 

- 0 ® © 


o 


& 




0 


Q 


O 


© 

© - © 

0 




lonx.iiAi'mcMAi’ Of maus 

2l\( tKUl 

Figure 2. Black spots show 61 "possible dust clouds" (defined in text) 
plus 27 other red clouds from selected sources. Triangles 
show "core areas" where major dust storms have formed. 
Dust storms appear associated with basins (dust reservoirs?) 
and borders of dark az*eas . 
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